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O’Shea,  Anthony  Thomas  (M.S.,  Mechanical  Engineering) 

Development  of  Terahertz  Rayleigh  Scattering  Diagnostic  for  a  Solid  Rocket  Exhaust  Plume 
Thesis  directed  by  Assistant  Professor  Andrew  D.  Ketsdever 

Characterization  of  rocket  plumes  has  been  difficult  due  to  the  lack  of  non-intmsive  real-time 
measurement  techniques.  This  research  proposes  a  new  technique  using  terahertz  radiation  and  Rayleigh 
scattering  off  of  particulates  in  the  plume.  Previous  techniques  have  involved  collecting  particles  on  a 
probe  or  mesh.  Others  have  used  Mie  scattering  using  lasers,  but  this  technique  is  limited  to  firing  the 
beam  straight  across  the  plume.  Terahertz  radiation  wavelengths  in  relation  to  alumina  particle  sizes  align 
with  the  Rayleigh  criterion,  which  requires  that  the  wavelength  of  light  be  at  least  ten  times  greater  than  the 
radius  of  the  particle.  Rayleigh  scattering  would  allow  for  changing  the  location  of  the  detector  and  even 
co-locating  it  with  the  source.  Unfortunately,  terahertz  sources  with  the  required  power  output  are 
currently  in  production.  Thus,  a  scaled  experiment  was  setup  using  microwaves  at  2.45  GElz  with 
aluminum,  brass,  and  stainless  steel  spheres  of  radii  ranging  from  4.76  mm  to  9.53  mm.  A  magnetron 
power  source  was  used  to  provide  a  power  range  of  200-2000  Watts  into  an  anechoic  chamber.  The 
forward  power  was  measured  and  then  passed  through  a  horn  antenna,  scattered  off  of  particles,  and 
collected  in  another  horn  antenna.  This  detected  power  was  measured  and  recorded.  The  same  test  was 
performed  for  no  spheres  inside  the  chamber  and  the  difference  was  considered  the  scattered  power.  These 
results  were  compared  with  predictions  obtained  from  a  Matlab  code  written  for  this  research.  From  this 
experiment,  it  was  shown  that  one  can  use  Rayleigh’s  equations  to  determine  particle  size  based  on 
scattered  power  detected  if  the  number  density  is  known.  However,  the  uncertainties  in  index  of  refraction 
and  scattering  efficiency  as  well  as  a  low  signal  to  noise  ratio  led  to  large  error  bars  and  discrepancies 
between  predictions  and  experimental  data.  Also,  a  notional  terahertz  diagnostic  instrument  was  developed 
that  would  operate  at  a  range  of  20  km  using  a  1  kW  source  and  a  3m  diameter  parabolic  reflector  antenna. 
Further  research  is  needed  to  characterize  the  properties  of  alumina  especially  the  refractive  index  and 
scattering  efficiency  to  ensure  accurate  measurements  from  the  terahertz  instrument. 
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CHAPTER  1 


INTRODUCTION 

This  study  seeks  to  further  the  development  of  a  non-intrusive  diagnostic  instrument  that  can  be 
used  to  characterize  alumina  particles  in  solid  rocket  plumes.  Particle  size  and  density  in  rocket  plumes  are 
difficult  to  measure  accurately.  Some  current  methods  involve  collecting  particles  on  a  screen  or  placing  a 
probe  into  the  plume  itself.1  While  these  methods  have  produced  results,  they  require  a  static  setup  and  can 
be  subject  to  contamination.  This  study  will  investigate  the  use  of  terahertz  radiation  and  the  principles  of 
Rayleigh  light  scattering  to  accurately  determine  particle  number  density,  size,  and  velocity.  Terahertz 
wavelengths  range  between  3000  jim  to  30  qm  and  the  particles  inside  the  rocket  plume  are  between  1  pm 
and  100  pm.  The  ratio  between  wavelength  and  particle  radius  satisfies  the  Rayleigh  light  scattering 
criterion  that  the  wavelength  must  be  at  least  ten  times  larger  than  the  particle  radius.  Therefore,  terahertz 
radiation  is  ideally  suited  to  be  used  in  the  development  of  this  diagnostic  tool. 

Since  the  dawn  of  the  rocket  age,  scientists  and  engineers  have  wanted  to  accurately  characterize 
the  performance  of  rockets  in  order  to  find  mechanisms  for  improvement.  Most  gains  in  performance  can 
be  realized  in  the  combustion  process  of  the  rocket  fuels  inside  the  combustion  chamber.  For  solid  rocket 
motors  in  particular,  diagnostics  of  properties  inside  the  combustion  chamber  are  especially  difficult 
because  the  intense  heat  and  pressure  would  destroy  most  measuring  tools.  Engineers  have  developed 
models  to  try  and  predict  the  chemical  processes  occurring  inside  the  combustion  chamber.  Flowever,  these 
models  need  high  fidelity  experimental  data  for  verification  and  validation.  The  solution  to  this  issue  may 
lie  in  the  understanding  of  the  rocket  plume  and  its  constituents.  If  one  can  accurately  measure  the  size, 
number  density,  velocity,  and  temperatures  of  the  particles  inside  the  solid  rocket  plume,  the  modeled 
combustion  process  can  be  adjusted  to  meet  targeted  plume  parameters  (i.e.  the  data  can  serve  as  an  end 
state  for  which  the  modeled  combustion  process  can  be  adjusted  to  meet). 

Measuring  these  characteristics  can  be  difficult  due  to  the  extreme  nature  of  the  plume 
environment.  Obtaining  real-time  measurements  adds  a  greater  challenge.  Light  scattering  has  the 
potential  to  provide  the  non-intrusive  measurement  tool  for  size,  number  density,  and  velocity  of  particles 
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in  real-time.  In  general,  particle  size  influences  the  efficiency  and  performance  of  the  rocket.  The  goal  is 
to  have  the  smallest,  lightest,  and  fastest  particles  generated  by  the  combustion  process,  leading  to  a  higher 
specific  impulse  and  thrust.  Different  propellant  compositions  produce  distinct  plume  particle 
characterizations.  By  adjusting  these  compositions,  based  on  measured  results,  performance  and  fuel  use 
can  be  optimized.  This  tool  could  be  used  to  measure  particle  sizes  and  therefore  efficiency  in  startup  and 
shut  down  of  solid  rocket  motors.  Also,  by  measuring  the  particle  size,  one  can  compare  the  performance 
of  different  combustion  processes. 

In  addition  to  the  possibility  of  improved  efficiency  through  analysis  of  the  plume  constituents, 
another  application  lies  in  the  identification  of  these  rockets  using  the  plume  signature.  With  a  detailed 
map  of  the  particles  generated  by  the  solid  rocket  motor,  the  manufacturer  and  origin  of  the  rocket  could  be 
identified.  Missile  defense  and  early  warning  systems  could  utilize  plume  signatures  to  potentially  identify 
any  enemy  threats. 

Light  scattering  offers  an  effective  real-time,  non-intrusive  diagnostic  tool  that  can  be  employed 
for  almost  any  solid  motor.  Limiting  the  use  of  light  scattering  to  Rayleigh  theory  offers  the  advantage  of 
sensor  and  source  mobility  due  to  the  dipole  structure  of  the  scattered  light.  This  advantage  would  also 
include  the  collocation  of  the  sensor  and  source.  Rayleigh  theory  stipulates  the  particle  size  must  be  ten 
times  smaller  than  the  wavelength  of  the  interrogating  light.  This  study  stems  from  the  fact  that  the  sizes  of 
particles  in  a  solid  rocket  exhaust  plume  and  the  wavelengths  of  terahertz  radiation  meet  this  criterion. 
Using  the  general  size  distribution  of  alumina  particles  that  occur  in  the  solid  rocket  motor  plume  coupled 
with  Rayleigh  theory,  a  notional  terahertz  diagnostic  instrument  was  developed  in  order  to  provide  a 
baseline  for  performance  that  an  instrument  would  need  to  be  successful  in  characterizing  the  solid 
constituents  of  the  plume.  Rayleigh  scattering  theory  was  combined  with  free  space  transmission  equations 
to  generate  a  single  equation  that  governs  the  ratio  between  transmitted  and  received  power.  From  this 
equation,  one  can  determine  if  a  real  system  would  have  sufficient  power  output  and  antenna  gain  to 
accurately  measure  the  scattered  return  from  particles  inside  the  solid  rocket  plume. 

In  order  to  test  the  feasibility  of  employing  Rayleigh  theory  to  the  diagnostic  development  for  the 
characterization  of  particles,  a  proof-of-concept  experiment  was  developed  involving  microwaves  at  2.45 
GHz  and  metal  spheres  of  several  sizes  and  materials.  The  experiment  was  scaled  based  on  the  diagnostic 
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wavelength  used  and  the  sphere  radius.  It  was  designed  to  analyze  the  possibility  of  using  Rayleigh 
scattering  to  measure  particle  size  while  determining  the  differences  between  theory  and  experiment.  It 
also  identified  some  of  the  complications  that  might  occur  when  using  a  terahertz  diagnostic  instrument. 

The  terahertz  region  of  the  electromagnetic  spectrum  between  100  GHz  and  10  THz  (wavelength 
range  from  3000  jim  to  30  fim)  occupies  a  large  portion  of  the  spectrum  between  the  infrared  and 
microwave  bands.  Compared  to  the  relatively  well-developed  science  and  technology  at  microwave, 
optical,  and  x-ray  frequencies,  developments  in  the  terahertz  band  are  limited  and  remain  relatively 
unexplored.2  The  shorter  wavelengths  at  terahertz  frequencies  allow  the  use  of  smaller  and  lighter 
components,  which  is  important  in  military  and  air-borne  applications.  In  addition,  compared  to  infrared 
and  optical  wavelengths,  atmospheric  attenuation  in  the  terahertz  region  is  relatively  low.3 

One  non-intrusive  light  scattering  technique  that  may  be  employed  with  the  particle  sizes  for  solid 
rocket  motors  is  Mie  scattering,  but  the  setup  would  have  to  be  static  due  to  the  general  forward-scattered 
direction  of  Mie  scattered  light.  Lasers  have  also  been  applied  to  hydrogen  and  oxygen  rocket  engines 
using  laser  and  Rayleigh  scattering.4  This  study  examines  solid  rocket  motors  which  produce  particles 
between  1  pm  and  100  pm.5  For  these  size  particles,  terahertz  wavelengths  from  10  pm  to  1  mm  are  ideally 
suited  to  measure  the  particles  in  the  plume  with  Rayleigh  scattering.  The  setup  is  not  restricted  to  a  static 
fire  and  could  be  used  in  a  real  time  launch  configuration.  The  use  of  Rayleigh  theory  also  removes  the 
large  complications  surrounding  the  application  of  Mie  theory  by  simplifying  the  intensity  ratio  to  a  basic 
sixth  order  polynomial  from  a  combination  of  several  associated  Legendre  polynomials  and  Riccati-Bessel 
functions. 

The  understanding  of  multiphase  rocket  exhaust  plume  flows  is  important  to  many  scientific  and 
engineering  applications.  Chemical  combustion  processes  involve  high  temperature,  multiphase  flows 
from  liquid  fuel  injectors  in  internal  combustion  engines  to  plume  effluents  at  the  exit  plane  of  a  rocket 
nozzle.  The  terahertz  region  of  the  electromagnetic  spectrum  may  have  a  large  impact  on  the  general 
understanding  of  these  two-phase  flows  for  particle  sizes  between  1  pm  and  100  pm.  Rayleigh  and  Raman 
spectroscopy  have  been  used  to  non-intmsively  diagnose  homogenous  and  heterogeneous  cluster 
formation,  condensation,  and  evaporation  in  a  high-temperature  gaseous  flow  environment.  Terahertz 
excitation  sources  bring  about  a  new  possibility  to  interrogate  flows  that  otherwise  could  not  be 
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investigated  with  normal  Rayleigh  or  Raman  techniques  due  to  the  correlation  between  terahertz 
frequencies  and  the  range  of  particle  size  with  respect  to  the  Rayleigh  criterion. 
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CHAPTER  2 


TERAHERTZ  SOURCES  AND  APPLICATIONS 
A.  Current  Sources:  Capabilities  and  Applications 

Terahertz  sources  have  seen  increased  development  in  the  past  decade  due  to  the  rising  interest  in 
the  wide  variety  of  applications  such  as  space,  security,  and  imaging.  As  mentioned  earlier,  the  terahertz 
regime  is  generally  held  to  be  between  100  GHz  and  10  THz.  This  is  a  region  where  high  power  sources 
have  been  difficult  to  develop  and  where  a  great  deal  of  research  is  currently  being  conducted.  Most 
terahertz  sources  today  range  in  output  power  from  single  micro-Watts  to  tens  of  Watts.6  7  There  are 
generally  two  ways  to  produce  terahertz  radiation:  lasers  and  gyrotrons.  Some  laser  handbooks  have 
compiled  the  output  powers  for  lasers  at  these  frequencies.8'9 10  From  this  data  and  using  the  general 
Rayleigh  theory  equation,  the  ratio  of  scattered  intensity  to  incident  intensity  was  determined  for  each  of 
the  different  laser  sources.  These  numbers  were  compared  with  the  predicted  performance  of  a  gyrotron 
source  that  will  be  discussed  later  in  this  chapter.  The  results  of  this  comparison  can  be  seen  below  in 
Figure  1. 
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Figure  1:  Intensity  ratio  for  several  sources 
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From  this  figure,  one  can  see  that  only  some  gas  lasers  and  free  electron  lasers  can  come  close  to  matching 
the  intensity  ratio  expected  from  the  gyrotron  source.  The  lasers  that  do  match  or  surpass  the  ratio  are  not 
considered  very  reliable,  and  it  will  become  apparent  that  a  high  output  power  will  be  required  for  a 
notional  terahertz  diagnostic  instrument  to  operate  effectively  for  this  application.  Also,  it  should  be  noted 
that  many  types  of  sources  cannot  attain  the  necessary  frequency  to  be  considered  in  the  terahertz  regime. 

While  current  terahertz  sources  may  not  have  large  output  powers,  many  applications  for  these 
relatively  low  power  sources  have  been  developed.  One  of  the  most  prominent  environments  for 
application  is  in  space.  Many  interstellar  particles  emit  radiation  in  the  terahertz  regime.  This  has  sparked 
great  interest  in  the  development  of  sensors  that  are  able  to  characterize  these  emissions.  This  information 
can  begin  to  answer  questions  concerning  the  composition  and  origins  of  the  universe.  Also,  terahertz 
sensors  can  help  characterize  the  atmospheric  properties  of  other  large  bodies  such  as  asteroids  and  comets. 
These  characteristics  can  lead  to  improvements  in  the  models  of  Earth’s  atmosphere.11  One  method  of 
employing  these  space  applications  would  be  with  ground  stations.  Gaidis12  discusses  the  spaced-based 
application  of  these  sensors  used  in  the  study  of  astronomy,  atmospheric  chemistry,  and  planetary  science. 

Another  large  area  of  terahertz  development  and  application  is  in  the  imaging  of  objects.  The  area 
has  a  broad  field  of  application  including  structural  analysis,  medical  procedures,  and  security  screening. 
Hu  and  Nuss13  demonstrated  the  ability  to  generate  terahertz  images  of  objects  such  as  semiconductors. 
They  used  terahertz  time -domain  spectroscopy  in  which  a  pulsed  laser  fired  into  a  terahertz  transmitter  that 
emitted  the  signal  across  the  test  specimen.  The  signal  was  collected  and  measured  by  a  terahertz  detector 
whose  data  was  resolved  into  an  image  of  the  test  piece.  Liu  et  al14  took  this  concept  and  applied  it  to 
defense  application  such  as  the  screening  of  packages  for  weapons  and  explosives.  They  demonstrated  this 
ability  by  creating  an  image  of  a  briefcase  with  a  knife  concealed  inside  and  the  structure  of  a  shoe. 
Bogue15  discussed  the  improvements  that  have  been  made  with  terahertz  applications.  He  explains  the  uses 
of  terahertz  radiation  in  the  medical,  security,  and  imaging  fields.  He  also  writes  about  the  new  security 
screening  for  airports  that  employs  terahertz  radiation.  This  study  seeks  to  add  another  application  for 
terahertz:  the  non-intrusive,  real-time  diagnostic  of  solid  rocket  exhaust  plumes. 
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B.  Gyrotrons:  Theory  and  Proposed  Source 

Gyrotron  devices  are  of  great  interest  because  of  their  ability  to  provide  high  average  and  high 
peak  power  at  millimeter  wavelengths  which  correspond  to  frequencies  from  30  to  300  GHz.  Current 
gyrotrons  have  increased  the  average  power  output  generated  by  radio  frequency  (RF)  sources  by  several 
orders  of  magnitudes.  In  general  gyrotrons  are  a  type  of  free  electron  maser.  A  gyrotron  consists  of 
vacuum  tubes  and  magnets  that  bunch  electrons  using  a  strong  magnetic  field.  A  cross-section  of  the 
electron  beam  for  a  gyrotron  can  be  seen  in  Figure  2. 


Basically,  some  electrons  are  decelerated  by  the  electric  field  which  increases  their  cyclotron 
frequency  while  others  are  accelerated  and  thus  their  cyclotron  frequency  decreases.  This  change  in 
cyclotron  frequencies  causes  the  electrons  to  bunch  causing  transverse  deceleration  from  which  energy  can 
be  extracted.  There  are  currently  several  high  power  gyrotron  systems;  one  of  which  is  a  Gyromonotron 
oscillator.  Gyrotron  oscillators  have  been  studied  since  the  late  1960s.  These  machines  are  capable  of 
long-duration  pulses  at  high  powers  which  effectively  place  them  in  a  continuous  wave  (CW)  regime.  One 
system  operates  at  170  GHz  with  an  output  power  between  0.5  MW  and  1  MW  for  durations  spanning  from 
0.5  to  multiple  seconds.  Another  state-of-the-art  device  is  the  gyrotron  driven  by  intense  relativistic 
electron  beams  (IREBs).  IREBs  which  have  very  high  current  and  high  beam  voltage  are  created  by  using 
explosive  emission  cathodes.  They  have  durations  less  than  nanoseconds  per  pulse.  A  gyro-traveling- 
wave-tube  (TWT)  mechanism  was  used  in  conjunction  with  IREBs  in  order  to  create  a  system  that  operates 
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at  35  GHz  with  an  output  power  of  100  MW.  Another  high  power  system  is  a  gyro-amplifier  which  can  be 
used  in  millimeter-wave  radar.  These  systems  amplify  the  signal  generated  by  a  gyro-TWT,  a 
gyroklystron,  or  a  gyro-twystron  to  obtain  high  power  at  certain  frequencies  for  radar  application.  Systems 
have  operated  around  95  GHz  with  output  power  ranging  from  2.5  kW  to  65  kW.  Gyrotons  have  also  been 
applied  to  electronic  counter  measures  (ECM).  In  these  systems  a  gyro-backward-wave-oscillator  (BWO) 
has  been  used.  In  ECM,  it  is  important  not  only  to  have  high-power  but  also  tunability.  Gyro-BWOs  offer 
this  capability  and  have  demonstrated  a  frequency  sweep  from  27-32  GHz  with  an  output  power  of  113 
kW.  Gyromonotron  oscillators  are  used  in  industry  and  have  shown  output  powers  around  200  kW  pulsed 
and  20  kW  CW  at  a  frequency  of  60  GHz.  Gyroklystrons  can  also  be  used  to  drive  accelerators.  Some 
systems  have  produced  output  powers  of  100  MW  at  19.76  GHz.16 

The  Electrical  and  Computer  Engineering  (ECE)  Department  at  the  University  of  Colorado  at 
Colorado  Springs  has  been  working  in  conjunction  with  this  study  to  develop  a  gyrotron  system  that  will 
operate  at  a  frequency  of  300  GHz  in  a  CW  regime  with  an  output  power  of  1  kW.  While  the  power  of  this 
proposed  source  does  not  match  the  state-of-the-art  systems  mentioned  above,  it  does  have  some  new 
characteristics  that  make  it  a  worthwhile  pursuit.  The  first  of  these  characteristics  is  the  operating 
frequencies.  300  GHz  lies  at  the  low-end  of  the  terahertz  regime.  This  presents  a  unique  challenge  because 
in  order  to  obtain  this  frequency  the  construction  process  will  have  to  involve  micro-fabrication  techniques. 
This  leads  to  the  second  major  advantage  of  this  system  which  is  the  compact  size.  Not  only  will  the 
system  itself  be  small,  but  the  necessary  magnet  will  also  be  substantially  smaller.  In  the  previous  systems, 
the  magnetic  field  required  is  rather  large,  thus  requiring  heavy  and  large  solenoids  to  generate  the 
appropriate  magnetic  field.  This  new  system  would  only  need  a  small  periodic  permanent  magnet.  The 
performance  characteristics  of  the  proposed  source  can  be  seen  below  in  Table  l.17 
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Device  Parameters 


Voltage 

40  kV 

Current 

2.0 -3.0  A 

a  =  v1/vz 

1.5 -2.0 

Avz/vz 

5% 

Magnetic  field 

1.1  Tesla 

Operating 

frequency 

>  300  GHz 

Output  power 

>  lkW 

Circuit  diameter 

~  700  pm 

Mode 

TEio,i(rr  mode) 

Harmonic  number 

s  =10 

Table  1:  Proposed  gyrotron  source  parameters17 


With  these  two  advances  in  gyrotron  technology,  the  new  source  would  be  beneficial  to  continuing  the 
work  of  this  study.  This  source  would  provide  the  basis  for  characterizing  the  Rayleigh  scattering  from 
alumina  particles  so  that  an  actual  terahertz  diagnostic  instrument  could  be  employed  for  the  real-time,  non- 
intmsive  diagnostics  of  a  solid  rocket  exhaust  plume. 
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CHAPTER  3 


SOLID  ROCKET  EXHAUST  PLUME  CHARACTERISTICS 

With  each  solid  rocket  there  are  different  chemical  make-ups  and  mixtures  which  lead  to  varying 
performance  in  terms  of  specific  impulse  and  total  thrust.  Most  solid  rocket  motors  use  three  basic 
ingredients  in  the  propellant  grain.  Those  ingredients  are  powdered  aluminum  as  a  fuel,  ammonium 
perchlorate  as  an  oxidizer,  and  hydroxyl-terminated  polybutadiene  as  a  binder  for  the  two  combustion 
components  with  varying  mixture  ratios.  The  differences  in  mixture  ratios  can  lead  to  a  variety  of  levels  of 
chemical  exhaust  products  and  affects  performance  of  the  rocket.  These  variances  in  the  rocket  properties 
need  to  be  understood  in  order  to  improve  performance  and  efficiency,  as  well  as  optimizing  prediction  and 
modeling  capabilities. 

A.  Rocket  Exhaust  Plume  Environment 

Inside  most  rocket  exhaust  plumes  there  exists  a  phenomenon  known  as  two-phase  flow.  After  the 
flow  has  expanded  and  cooled,  gas  molecules  as  well  as  solid  particles  constitute  the  flow  of  the  plume. 
The  ratio  of  solid  particles  to  gas  varies  with  every  rocket  depending  on  the  propellant  and  mixtures.  Also, 
the  effect  on  the  flow  properties  stems  back  to  the  varying  ratios.  In  liquid  hydrocarbon  fuels,  the  solid 
phase  is  composed  of  carbon  soot  particles.  In  solid  rocket  motors  using  aluminum  powder,  the  solid  phase 
is  comprised  of  aluminum  oxide.  The  aluminum  oxide  particles  are  generally  larger  than  soot  particles.  A 
comparison  of  the  two  rocket  types  and  a  general  description  of  the  two  phase  flow  environment  can  be 
seen  in  Figure  3. 18 
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Figure  3:  The  two-phase  flow  environment  for  liquid  and  solid  rockets18 

In  order  to  characterize  the  plume  environment,  several  prediction  codes  have  been  developed. 
One  code  is  the  SPF/SIRRM  which  was  designed  to  model  the  “nozzle  flow-field,  body  flow-field, 
base/separated  region,  and  the  plume  flow-field.”  This  code  also  models  the  radiation  intensity  generated 
by  a  rocket  which  has  missile  defense  applications  for  detection  and  identification.  However,  deficiencies 
with  the  code  lead  to  a  poor  handling  of  initial  conditions  in  the  nose  region  of  the  rocket.  Another  code  is 
CHARM  which  is  a  compilation  of  several  codes.  The  CHARM  code  performs  similar  functions  as  the 
SPF/SIRRM  code  by  modeling  the  particles  in  the  flow  of  a  rocket  and  determining  the  predicted  emissive 
properties  of  plumes.  One  shortcoming  of  the  CHARM  code  is  that  it  has  difficulty  accurately  modeling 
some  reaction  rates  in  the  shock  layer.19  Numerical  computations  were  also  performed  on  a  two-phase  flow 
by  Owis  and  Hashem.20  They  analyzed  the  effects  of  two-phase  flow  in  satellite  control  thrusters.  They 
numerically  solved  the  unsteady  compressible  Navier-Stokes  equation  to  determine  the  effects  of 
condensing  gaseous  nitrogen.  They  found  that  in  general  this  condensation  lowered  the  thrust  of  the 
satellite  control  thruster.  If  a  terahertz  diagnostic  instrument  were  built  and  implemented,  some  of  the 
unknown  parameters  used  in  these  codes  could  be  better  quantified  leading  to  higher  fidelity  predictions  of 
particle  size  and  number  density  inside  the  exhaust  plume. 

B.  Aluminum  Oxide  Properties 

Aluminum  oxide  is  the  primary  constituent  of  concern  for  this  study.  After  the  flow  has  expanded 
and  cooled  inside  the  solid  rocket  exhaust  plume,  aluminum  oxide  exists  in  the  form  of  solid  particles  and 
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is  sometimes  referred  to  as  alumina  or  synthetic  sapphire.  If  the  notional  terahertz  instrument  is  ever  to  be 
used  in  a  real-time  diagnostic  of  a  solid  rocket  exhaust  plume,  many  characteristics  of  alumina  will  have  to 
be  taken  into  account.  The  first  of  those  characteristics  will  be  the  refractive  index  of  alumina  at  the 
specific  operating  frequency.  Krishnan  el  al 21  performed  an  experiment  to  find  the  refractive  index  of 
liquid  aluminum  oxide.  They  used  a  process  called  laser  ellipsometry  which  measures  laser  light  reflected 
off  a  sample  and  creates  a  temperature  difference  between  the  two  sides  of  the  sample.  Their  experiment 
found  only  the  real  part  of  the  refractive  index.  The  wavelength  was  specific  at  0.6238  pm  which  yielded  a 
refractive  index  of  1.744  ±  0.016.  This  value  was  used  to  design  the  notional  terahertz  diagnostic 
instmment.  It  should  be  noted  that  the  complex  index  of  refraction  is  of  important  value  as  well  and  is  used 
in  Mie  theory  calculations.  Gupta  and  Wall22  performed  an  experiment  using  a  0.5  pm  wavelength  and 
found  both  the  real  and  imaginary  parts  of  the  refractive  index.  They  reported  a  complex  value  of 
1.77+0.005/  which  illustrates  the  presence  of  the  imaginary  component.  It  should  also  be  taken  into 
account  that  there  is  an  apparent  dependence  on  wavelength  based  on  the  previous  two  reported  indices. 
Gryvnak  and  Burch23  conducted  experiments  to  try  to  relate  the  index  of  refraction  to  both  wavelength  and 
temperature.  Malitson,  Murphy,  and  Rodney24  also  looked  at  the  same  relation,  but  at  lower  temperatures. 
They  generated  the  following  data  in  Table  2  for  various  wavelengths  and  temperatures. 


Wavelength  (pm) 

17.0°C 

24.0°C 

31.0°C 

0.7065 

1.76294 

1.76303 

1.76312 

0.6438 

1.76538 

1.76547 

1.76556 

0.5791 

1.76862 

1.76871 

1.76880 

0.5770 

1.76814 

1.76883 

1.76892 

0.5461 

1.77068 

1.77077 

1.77086 

0.4358 

1.78110 

1.78120 

1.78130 

0.4047 

1.78572 

1.78582 

1.78592 

Table  2:  Index  of  refraction  of  synthetic  sapphire  at  visible  wavelengths  and  three  temperatures24 

Using  this  data  combined  with  the  results  from  Gryvnak  and  Burch,  Plass25  was  able  to  generate  data  for 
refractive  index  of  aluminum  oxide  at  higher  temperatures.  The  real  part  of  the  index  was  assumed  to  be 
increased  by  2.9xlOA-5/°C.  The  results  are  shown  in  Table  3. 
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Wavelength 

Til 

Temperature 

M 

1200°C 

1500°C 

1600°C 

1700°C 

2020’C 

0.5 

1.805 

1.814 

1.817 

1.82 

1.83 

1 

1.785 

1.794 

1.797 

1.80 

1.81 

2 

1.775 

1.784 

1.787 

1.79 

1.80 

3 

1.745 

1.754 

1.757 

1.76 

1.77 

4 

1.715 

1.724 

1.727 

1.73 

1.74 

5 

1.665 

1.674 

1.677 

1.68 

1.69 

6 

Wavelength 

1.575 

1.584 

1.587 

m 

Temperature 

1.59 

1.60 

M 

120O“C 

1500°C 

1600 “C 

1700°C 

2020°C 

0.5 

1.8  X  10-' 

5.2  X  10“ 

1.1  X  10“ 

1.6  X  10“ 

1.2  X  10“ 

1 

6.0  X  10-* 

4.4  X  10“ 

1.1  X  10“ 

1.8  X  10“ 

1.6  X  10“ 

2 

5.3  X  10“ 

3.8  X  10“ 

1.0  X  10“ 

1.9  X  10“ 

2.2  X  10“ 

3 

1.9  X  10-’ 

8.4  X  10“ 

1.7  X  10“ 

2.4  X  10“ 

3.6  X  10“ 

4 

1.4  X  10~‘ 

2.0  X  10“ 

2.4  X  10“ 

2.9  X  10“ 

3.5  X  10“ 

5 

1.6  X  10-‘ 

2.3  X  10“ 

2.8  X  10“ 

3.2  X  10“ 

3.7  X  10“ 

6 

9.1  X  10“ 

1.1  X  10“ 

1.4  X  10“ 

1.6  X  10“ 

1.9  X  10“ 

Table  3:  Index  of  refraction  dependence  on  wavelength  and  temperature25 


In  the  above  table  nj  and  n2  represent  the  real  and  imaginary  components  of  the  complex  index  of 
refraction.  From  these  two  tables,  it  is  obvious  that  before  any  diagnostic  of  solid  rocket  exhaust  plumes 
can  be  performed  using  the  notional  terahertz  diagnostic  instrument,  the  optical  properties  of  alumina  in  the 
terahertz  regime  must  be  fully  characterized.  The  refractive  index  must  be  determined  for  the  operating 
wavelengths  as  well  as  the  temperature  dependence  characterized  for  terahertz  radiation.  The  refractive 
index  plays  an  important  part  in  the  determination  of  the  scattering  cross-section  which  is  part  of  the 
determination  of  the  power  ratio  in  both  Mie  and  Rayleigh  theories.  No  known  values  for  the  index  of 
refraction  of  alumina  at  terahertz  wavelengths  and  temperatures  inside  the  plume  have  been  found  in  the 
literature.  These  values  would  be  needed  for  the  successful  operation  of  a  terahertz  diagnostic  instalment. 
C.  Measurements  and  Predictions  of  Aluminum  Oxide 

Over  the  years,  many  solid  rocket  launches  have  yielded  data  on  the  particles  generated  from  the 
combustion  process.  One  prominent  paper  collected  data  from  many  sources  in  order  to  compile  data  about 
aluminum  oxide  particle  size.  R.W.  Hermsen26  used  data  generated  between  1962-1978  to  construct  a  large 
table  of  data  including  60  different  motor,  propellant,  and  chamber  pressure  varieties.  The  complete  set  of 
data  from  this  comparison  can  be  seen  in  Appendix  A.  From  this  table  along  with  other  prediction  data, 
Flermsen  was  able  to  generate  a  single  chart  which  relates  the  mean  aluminum  oxide  particle  diameter  to 
the  throat  diameter  of  the  solid  rocket  nozzle  as  in  Figure  4. 
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Throat  diameter  (in.) 

Figure  4:  Comparison  of  particle  size  correlations26 

One  type  of  aluminum  oxide  particle  measurement  technique  calls  for  inserting  a  probe  into  the 
exhaust  plume  for  a  known  amount  of  time  and  then  measuring  the  collected  particles.  Laredo  et  al11 
collected  particles  in  two  different  ways.  One  was  scraping  the  particles  from  the  walls  of  the  nozzle  and 
combustion  chamber  while  the  other  involved  a  stainless  steel  wedge  covered  in  a  copper  sheet.  This  probe 
used  a  flat  front  with  tapered  sides  and  was  inserted  into  the  exhaust  plume  for  0.5s.  Once  the  collection 
was  completed,  the  sample  was  immediately  analyzed  using  a  scanning  electron  microscope.  Afterwards, 
the  data  was  complied  and  graphed. 

Other  engineers  have  strived  to  generate  new  measurements  and  measurement  techniques  to 
determine  aluminum  oxide  particle  size.  J.  K.  Sambamurthil  developed  a  new  probe  device  to  enter  into 
the  exhaust  plume  environment  and  collect  the  particles  as  the  rocket  is  firing  in  a  static  position.  He 
developed  a  system  that  launched  collection  darts  through  the  exhaust  plume.  Each  of  the  darts  was 
covered  with  a  copper  tape  that  collected  the  particles  during  its  duration  in  the  plume.  After  the  rocket  had 
been  fired  and  the  data  samples  collected,  the  copper  tapes  from  the  dart  were  analyzed  using  a  scanning 
electron  microscope.  Figure  5  shows  images  of  a  typical  collection  tape  when  viewed  under  a  microscope. 
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Figure  5:  Typical  collection  images  from  scanning  electron  microscope] 


From  the  data,  Sambamurthi  concluded  that  the  particles  collected  are  mostly  spherical  and  have  a  variety 
of  different  sizes.  This  dart  system  and  its  data  may  prove  to  be  beneficial  in  the  validation  stages  of  an 
actual  terahertz  diagnostic  system. 

in  addition  to  aluminum  oxide  particle  measurements,  the  prediction  of  the  particles’ 
characteristics  based  on  the  properties  of  the  solid  rocket  motor  has  become  a  field  of  interest.  O.  B. 
Kovalev"8  has  developed  prediction  methods  for  motor  and  plume  particle  size.  This  prediction  model 
relies  on  some  previous  models,  but  in  general  creates  a  new  model  which  advances  the  modeling  of 
agglomerate  particles  that  enter  the  flow.  He  validates  his  model  against  the  data  generated  by  R.W 
Hermsen’s  study  and  finds  that  his  model  is  in  line  with  historical  data.  Some  other  models  for  particle  size 
distribution  and  particle  density  include  the  model  developed  in  Reference  29  where  the  author  develops  a 
model  based  on  historical  data  collected  from  actual  launches  at  altitudes  in  the  stratosphere.  Using  this 
data,  he  adjusted  the  model  to  generate  similar  results  for  particle  size  and  number  density.  These  models 
could  be  validated  by  data  obtained  by  an  actual  terahertz  diagnostic  system. 


Distribution  A:  Approved  for  public  release.  Distribution  unlimited. 


CHAPTER  4 


PREVIOUS  RAYLEIGH  EXPERIMENTS 

In  theory,  Rayleigh  scattering  experiments  can  be  performed  with  any  kind  of  electromagnetic 
radiation  or  light.  The  only  consideration  that  must  be  made  when  trying  to  perform  a  Rayleigh  scattering 
experiment  is  the  size  of  the  particles.  This  property  of  the  particles  will  dictate  what  wavelengths  can  be 
used  for  Rayleigh  scattering.  The  criterion  is  that  the  particle  radius  must  be  at  least  ten  times  smaller  than 
the  wavelength  of  the  interrogating  light. 

A.  General  Rayleigh  Experiments 

In  the  past,  there  have  been  many  experiments  to  measure  very  small  molecules  and  clusters  of 
atoms.  Since  the  small  size  allows  for  a  small  wavelength  while  still  satisfying  the  Rayleigh  criterion,  there 
have  been  several  experiments  using  lasers  which  tend  to  have  wavelengths  in  the  100  to  1000  nanometer 
range.  In  Reference  30,  the  investigators  were  able  to  demonstrate  the  capability  of  measuring  the  size  of 
argon  clusters  using  Rayleigh  scattering.  In  the  experiment,  they  created  atomic  clusters  by  flowing  argon 
gas  through  a  nozzle  into  a  vacuum  environment  using  adiabatic  expansion.  In  order  to  further  simplify  the 
Rayleigh  scattering  theory,  they  related  the  number  density  of  the  argon  clusters  to  the  reservoir  pressure 
assuming  that  all  of  the  argon  had  formed  into  clusters.  The  experiment  used  a  pulsed  Nd:YAG  laser  at  a 
wavelength  of  532  nm  to  generate  the  Rayleigh  scattered  light.  The  laser  output  was  less  than  1  mJ  with  a 
15  ns  pulse.  The  stray  reflected  light  was  absorbed  in  order  to  prevent  interference  with  the  measurements. 
The  results  from  the  experiment  can  be  seen  in  Figure  6. 
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Figure  6:  Dependence  of  argon  cluster  size  on  the  backing  pressure30 


The  authors  were  successful  in  showing  through  Rayleigh  scattering  that  the  size  of  atomic  argon  clusters  is 
directly  dependent  upon  the  reservoir  pressure. 

Seasholtz,  Zupanc,  and  Schneider31  also  use  Rayleigh  scattering  as  a  diagnostic  tool.  They  show 
the  capability  to  measure  gas  density,  temperature,  and  velocity  in  a  hydrogen-oxygen  rocket  exhaust 
plume.  Their  diagnostic  focuses  on  resolving  the  Rayleigh  scatter  spectrum.  In  order  to  accomplish  this 
task,  they  utilize  a  Fabry-Perot  interferometer.  The  experiment  was  performed  using  a  100  N  thruster  with 
a  chamber  pressure  of  690  kPa  and  an  exit  pressure  of  2.1  kPa.  These  pressures  generated  a  mean  free  path 
of  about  20  pm.  This  led  the  team  to  apply  collisionless  Rayleigh  scattering  theory.  The  source  used  was  a 
1  W  argon-ion  laser  with  a  wavelength  of  514.5  nm.  The  laser  was  fired  into  the  rocket  exhaust  plume  and 
then  terminated  after  the  plume.  The  scattered  light  was  collected  and  sent  through  the  Fabry-Perot 
interferometer  and  then  into  a  photo-multiplier  tube.  The  tests  ran  for  100s  after  which  the  collected  data  is 
processed  and  then  plotted.  Figure  7  shows  the  results  of  the  experiments  with  different  oxidizer  to  fuel 
ratios  compared  with  JANNAF-TDK  code. 
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Figure  7:  Results  from  spectrally  resolved  Rayleigh  scattering  experiment31 


Seasholtz,  Zupanc,  and  Schneider  also  discuss  some  sources  of  error  that  arose  in  their  experiment.  One  of 
the  most  prominent  sources  were  thermal  broadening  coupled  with  spectral  broadening  of  the  Rayleigh 
scattered  light  which  led  to  an  increase  in  the  measured  temperature.  Another  error  source  was  excess 
hydrogen  gas  in  the  flow  which  increased  the  measured  number  density.  Such  factors  will  also  have  to  be 
addressed  in  the  application  of  terahertz  diagnostics  to  solid  rocket  exhaust  plumes. 

Zupanc32  continued  research  in  this  area  and  performed  a  follow-on  experiment  that  involved 
using  pulsed  laser  light  to  interrogate  a  plume.  He  focused  on  determining  the  exit  plane  velocity  for  a 
small  scale  hydrogen-oxygen  rocket.  He  relied  on  Rayleigh  scattering  light  and  the  Doppler  shift  of  the 
scattered  light  in  order  to  make  his  measurements.  The  data  was  compared  with  the  predictions 
RK/RPLUS  code  which  is  a  full  Navier-Stokes  model.  The  radial  results  were  in  good  agreement,  but  there 
were  large  errors  associated  with  the  axial  results. 
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Seasholtz  also  continued  his  use  of  Rayleigh  scattering  as  a  diagnostic  technique  along  with 
Buggele  and  Reeder33.  The  new  experiment  was  designed  to  again  use  Rayleigh  scattered  light  along  with 
a  Fabry-Perot  interferometer  to  perform  flow  diagnostics  for  a  supersonic  wind  tunnel.  The  main 
motivation  for  the  experiment  lies  in  the  non-intrusive  nature  of  a  Rayleigh  scattering  diagnostic. 
Conventional  probes  affected  the  flow  too  much  to  produce  accurate  results.  Therefore,  the  team  used  a 
Nd:YAG  laser  at  a  wavelength  of  532  nm.  The  laser  was  fired  into  the  flow  in  single  bursts  and  collected 
at  a  scattering  angle  of  88.9  degrees.  The  data  was  then  post-processed  to  determine  the  flow 
characteristics.  The  team  assumed  that  the  flow  was  adiabatic  for  the  generation  of  their  results  which  can 
be  seen  in  Table  4.  The  measurements  were  made  at  two  separate  points  in  the  flow:  in  front  of  and  behind 
a  shock. 


Parameter 

RDG  1298 

RDG  1300 

Mass  flow  rate 

13  lb  ms"1 

16  lb  ms"1 

Total  temperature  T7  (plenum) 

296K 

296K 

Total  pressure  P7  (plenum) 

16.9  psia 

21.5  psia 

Static  pressure  P8 

6.28  psia 

7.37  psia 

Mach  number  MW  (wedge  probe) 

1.95 

1.95 

Mach  number  MN3 

0.94 

1.00 

Mean  velocity  (std  dev)  -  Rayleigh 

414  (51)  ms"1 

445  (42)  ms"1 

Mean  temperature  (std  dev)  -  Rayleigh 

209  (21)K 

196  (18)  K 

Mean  Mach  number  (std  dev)  -  Rayleigh 

1.44  (0.26) 

1.60  (0.22) 

Table  4:  Results  from  Rayleigh  scattering  flow  diagnostics33 

One  problem  that  was  encountered  was  the  long  amount  of  time  that  was  required  to  post-process  the  data. 
The  authors  concluded  that  a  faster  method  is  needed  for  this  diagnostic  technique. 

Seasholtz  in  conjunction  with  Panda34  extended  this  work  to  characterize  the  velocity  and 
temperature  inside  a  supersonic  free  jet  using  spectrally  resolved  Rayleigh  scattered  light.  The  approach  in 
the  experiment  utilizes  a  CW  laser  along  with  a  Fabry-Perot  interferometer  and  photo  multiplier  tubes.  The 
power  source  was  a  5  W  Nd:YV04  laser  that  operated  at  532  nm.  The  light  was  transmitted  into  the  jet 
stream  and  then  collected  and  sent  through  the  interferometer  and  the  photo  multiplier  tubes.  With  this 
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setup,  the  team  was  able  to  demonstrate  the  capability  to  measure  temperature  and  velocity  in  both  axial 
and  radial  direction  inside  a  Mach  0.99  and  Mach  1.395  jet  stream. 

Another  experiment  that  employed  Rayleigh  scattering  to  determine  the  size  of  a  substance  was 
conducted  by  Zhang  et  a/.35  In  this  test,  the  researchers  set  out  to  determine  the  size  of  glycogen  that 
agglomerates  in  a  solution.  They  use  a  laser  spectrometer  that  uses  a  200  mW  argon-ion  laser  at  514.5  nm. 
A  spectrofluorometer  was  also  used  to  measure  the  light  scattering  intensity  of  the  solution.  Using  these 
two  instruments,  the  team  was  able  to  show  the  presence  of  glycogen  in  a  solution  and  measure  the  average 
hydro-dynamic  radius  of  the  glycogen.  These  Rayleigh  diagnostic  experiments  lend  credibility  to  the 
feasibility  of  successfully  conducting  the  proof-of-concept  experiment  and  the  application  of  the  notional 
terahertz  diagnostic  instrument  to  the  characterization  of  solid  rocket  exhaust  plumes. 

B.  Microwave  Rayleigh  Experiments 

Previous  experiments  involving  microwave  frequencies  and  using  the  principles  of  Rayleigh 
scattering  are  of  great  interest  to  this  study  given  the  similar  nature  of  the  proof-of-concept  experiment. 
Many  of  these  experiments  involve  a  diagnostic  of  a  plasma  which  while  different  from  strictly  particles, 
still  provides  insight  into  the  workings  and  complications  surrounding  the  use  of  Rayleigh  scattering. 
Shneider  and  Miles36  began  using  Rayleigh  scattering  with  microwave  frequencies  to  investigate  the 
properties  of  small  plasma  objects.  Their  study  developed  a  method  that  could  be  used  as  a  diagnostic  for 
small  scale  plasmas  such  as  laser  sparks,  avalanche-streamer  transitions,  and  resonance-enhanced  multi¬ 
photon  ionizations  processes.  They  treated  a  plasma  as  a  source  of  scattered  dipole  radiation  when  it 
encounters  incident  microwave  radiation.  This  principle  sets  up  the  use  of  Rayleigh  scattering  theory  to 
predict  the  intensity  of  the  scattered  light  and  the  scattering  cross-section  of  the  plasma.  Once  they 
developed  their  equations  for  the  total  averaged  power  radiated  by  the  dipole  and  the  scattering  cross- 
section,  they  generated  plots  of  the  predicted  cross-section  of  the  plasma  using  simulated  incident  intensity 
at  frequencies  of  3  and  6  GHz.  The  predicted  results  can  be  seen  in  Figure  8. 
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The  methods  developed  can  be  used  to  measure  the  electron  number  density  and  electron-loss  rate  in  the 
plasma  based  on  the  Rayleigh  scattered  signal. 

In  Reference  37,  the  authors  apply  the  theories  developed  in  Reference  36.  Their  study  worked  to 
conduct  diagnostics  of  laser-induced  avalanche  ionization  in  air  using  microwaves  and  Rayleigh  scattering. 
The  team  first  develops  equations  for  predicting  the  results  from  the  experiment.  They  develop  formulas 
for  the  scattering  cross-section  of  the  laser  spark.  Once  the  formulas  are  obtained,  graphs  are  generated  for 
the  predicted  number  of  electrons  in  the  laser  spark.  Once  predictions  are  finalized,  the  experiment  is  set 
up.  In  order  generate  the  laser  breakdown  in  air,  a  Nd:YAG  laser  was  used.  For  the  scattering  setup  a 
tunable  Gunn-diode  power  source  was  used  with  an  output  power  of  15  mW  at  a  frequency  of  12.6  GHz. 
The  power  was  transmitted  through  a  WR75  horn  which  was  30  cm  away  from  the  laser  breakdown.  After 
the  light  had  scattered  it  was  collected  in  three  WR75  receiving  horns.  The  recorded  data  was  processed 
and  was  found  to  be  in  good  agreement  with  the  predicted  results.  Thus,  microwaves  and  Rayleigh 
scattering  were  shown  to  be  a  useful  diagnostic  technique  for  laser-induced  avalanche  ionization  in  air. 

In  Reference  38,  the  investigators  continued  their  research  using  microwaves  and  Rayleigh 
scattering  in  order  to  detect  resonance-enhanced  multiphoton  ionization  in  argon.  Their  first  steps  were  to 
generate  the  theory  for  predicting  the  amount  of  microwave  scattering  expected  from  the  ionization  in 
argon.  Once  the  theory  had  been  established,  they  set  up  an  experiment  use  a  Ti:sapphire  laser  at  261.27 
nm  to  create  the  ionization.  To  generate  and  measure  the  scattering,  a  10  mW  microwave  source  at  12.6 
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GHz,  which  corresponds  to  a  wavelength  of  2.3  cm,  was  used  along  with  WR75  horns.  After  the  test  was 
completed  the  recorded  data  was  compared  with  the  predicted  results  which  can  be  seen  in  Figure  9. 


Figure  9:  Experimental  and  predicted  results  of  microwave  Rayleigh  scattering  experiment38 


One  can  see  from  Figure  9  that  the  experimental  results  and  theory  are  in  good  agreement.  Thus,  the  team 
was  able  to  conclude  that  with  their  diagnostic  technique  they  can  successfully  detect  resonance-enhanced 
multiphoton  ionization  in  argon  using  Rayleigh  scattering. 
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CHAPTER  5 


THEORY  OF  APPLIED  LIGHT  SCATTERING 
A.  A  Comparison  of  Mie  and  Rayleigh  Theories 

1.  Mie  Scattering  Theory 

Mie  theory  is  an  elastic  light  scattering  theory  that  was  developed  by  Gustav  Mie  and  has 
sometimes  been  referred  to  as  Mie-Lorenz  theory.  Mie  started  with  Maxwell’s  equations  and  resolved 
them  into  equations  that  related  the  scattered  light  to  incident  light  as  well  as  defined  cross-sections  for 
scattering  and  absorption.  The  equations  that  are  most  pertinent  to  this  study  involve  the  relation  between 
the  scattered  intensity  from  a  spherical  scattering  center  to  the  incident  light  that  interacts  with  the 
scattering  center.  Van  de  Hulst39  has  given  a  concise  summary  of  this  theory,  and  the  final  equation  for  a 
single  spherical  scattering  particle  can  be  seen  below 

4(4  +  4) 


4  = 


2  k2d2 


k  is  related  to  the  wavelength  of  the  incident  light  by 

2n 

k  =  T 

The  two  i  terms  are  defined  by  the  amplitude  functions  for  the  scattered  light  by 
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These  amplitude  functions  can  be  defined  as 
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The  coefficients  a  and  b  are  defined  using  Riccati-Bessel  functions  by  the  equations 

<p'n(y)<Pn(x)  -  m(pn{y)(p'n{x) 
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m(p'n (y)(pn (x)  -  (pn C y)(p'n (x) 
b  — - 

n  m(p'n  (y)Cn (x)  -  <pn (y)C CO 

The  functions  cp  and  £  are  defined  by  Bessel  functions  by  the  following 


(pn(z)  = 


CnO)  = 


7TZ 


(2) 


2  ^n  +  l/2 


(z) 
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(9) 

(10) 


The  values  for  x  and  y  are 


x  =  kr 

(11) 

y  —  mkr 

(12) 

The  coefficients  n  and  r  from  Eq.  (5)  and  Eq.  (6)  are  represented  by  associated  Legendre  polynomials  by 
the  following  relations 


1  , 

nn(cos9)  = - -Pj}(cos9) 

sin0 

d  i 

f„(cos  9)  =  —  PJ  (cos  9) 


(13) 

(14) 


Using  this  theory  for  the  relation  between  incident  light  and  scattered  light,  one  can  determine  the  size  of  an 
arbitrary  spherical  particle  regardless  of  the  wavelength  of  the  incident  light. 

2.  Rayleigh  Scattering  Theory 

Rayleigh  scattering  is  attributed  to  Sir  John  Rayleigh  who  is  well  known  for  describing  the  blue 
color  of  the  sky.  He  concluded  that  the  ratio  of  scattered  light  intensity  to  incident  light  intensity  is 
proportional  to  the  sixth  power  of  the  particle  radius  and  inversely  proportional  to  the  distance  from  the 
particle  squared  and  the  fourth  power  of  the  wavelength  of  light.  Rayleigh  scattering  is  elastic  scattering 
where  the  frequency  of  scattered  light  is  equivalent  to  the  frequency  of  the  incident  light.  Fundamentally, 
Rayleigh  scattering  is  an  approximation  of  the  overarching  theory  of  Mie  scattering  which  is  also  elastic  in 
nature.  However,  there  is  one  important  limitation  on  Rayleigh  scattering.  The  particle’s  radius  must  be  at 
least  ten  times  smaller  than  the  wavelength  of  the  incident  light.  If  the  particle  is  any  larger  the  theory  does 
not  conform  to  the  values  obtained  using  Mie  theory.  An  advantage  of  Rayleigh  scattering  over  Mie 
scattering  is  the  relative  simplicity  in  calculating  predicted  scattered  powers.  Mie  theory  is  a  complicated 
process  involving  either  several  infinite  sums  of  Legendre  polynomials  or  multiple  infinite  sums  of  Bessel 
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functions  of  several  kinds  and  Hankel  functions.  Rayleigh  theory  simplifies  the  Mie  theory  equations  into 

one  simple  ratio  of  scattered  intensity  (7S)  to  incident  intensity  (/0)40'41 

2 


167r4r6 


(mz  —  IN 
\m2  +  2 ) 
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The  scattered  intensity  represented  in  Eq.  (15)  is  the  total  intensity  of  light  scattered  from  the  particle.  If 
one  wants  to  detect  scattered  light,  only  a  portion  of  this  light  will  end  up  at  a  given  sensor.  The  solid  angle 
of  the  sensor’s  collection  aperture  must  be  taken  into  account  using  the  following  equation. 
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This  integral  can  be  used  to  determine  either  the  mean  size  of  particles  or  the  number  density  of  particles 
with  a  certain  size. 

Another  property  of  light  scattering  in  the  Rayleigh  regime  is  the  scattered  intensity  has  an  almost 
uniform  dispersion.  In  the  Mie  regime  the  scattered  intensity  tends  to  be  in  the  same  direction  as  the 
incident  light.  This  fact  gives  an  experiment  using  the  Rayleigh  regime  of  scattered  light  more  flexibility  in 
set  up  because  more  light  can  be  easily  detected  in  a  wide  variety  of  sensor  locations.  Figure  10  illustrates 
the  difference  in  scattering  patterns  between  Mie  and  Rayleigh. 


Rayleigh  Scattering  Mie  Scattering  Mie  Scattering  for  larger  particles 


4*  Direction  of  Incident  Light 


Figure  10:  Depiction  of  elastic  light  scattering 


B.  Antenna  Theory 

An  important  characteristic  of  any  transmitting  device  is  the  antenna  used  to  transmit  the  radiation. 
An  antenna  offers  a  method  of  directionalizing  the  emitted  radiation  in  order  to  combat  the  large  losses  that 
are  accrued  when  propagating  through  free  space.  This  directionality  is  defined  in  terms  of  an  antenna’s 
gain.  A  high  gain  corresponds  to  a  very  acute  directional  ability,  while  a  small  gain  tends  toward  limited 
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directional  capabilities.  There  are  many  different  antenna  types  that  range  from  simple  dipole  antennas 
with  very  small  gain  and  therefore  very  limited  directional  ability  to  massive  parabolic  reflectors  such  as 
the  deep  space  communications  antennas  which  offer  very  high  gain  and  are  therefore  very  directional  in 
nature.  Each  different  antenna  has  its  advantages  and  disadvantages  for  its  specific  purpose. 

In  this  study  to  develop  a  diagnostic  tool  for  using  terahertz  radiation  to  interrogate  a  solid  rocket 
exhaust  plume,  two  different  antenna  types  were  investigated.  For  the  notional  terahertz  diagnostic 
instrument,  a  parabolic  reflector  was  selected  due  to  the  ability  to  achieve  the  high  gains  required.  For  the 
proof-of-concept  experiment  using  microwaves,  horn  antennas  were  selected  for  the  ease  of  use  and 
procurement.  Both  antenna  types  will  be  discussed  below. 

1.  Light  Propagation  Considerations 

Anytime  light  is  generated  or  emitted  it  must  travel  through  free  space.  The  intensity  of  the  light, 
and  thus  its  power,  is  inversely  proportional  to  the  square  of  the  distance  traveled.  This  natural  occurrence 
of  loss  is  deemed  free-space  loss  and  is  an  issue  that  must  be  dealt  with  in  every  kind  of  transmitting  piece 
of  equipment  from  radios  to  televisions  and  especially  satellite  communications.  To  combat  these  losses, 
engineers  and  scientists  have  developed  antennas  to  create  a  directional  signal  and  allow  more  of  the 
original  transmitted  power  to  be  preserved.  However,  in  order  to  know  what  antenna  to  use  and  how 
directional  it  must  be,  the  system  losses  must  be  understood.  The  free  space  loss  was  characterized  using 
the  following  equation.42 


This  would  account  for  the  free  space  loss  of  light,  but  another  consideration  must  also  be  acknowledged  in 
terms  if  attenuation  of  a  signal:  atmospheric  interference.  Many  studies  have  looked  at  the  problem  of 
atmospheric  attenuation  and  developed  graphs  depicting  the  atmospheric  windows  for  transmitted  radiation 
while  minimizing  losses.43  Figure  11  below  shows  one  graph  that  highlights  the  specific  region  in  which 
the  terahertz  source  in  development  would  operate.44 
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Figure  11:  Atmospheric  attenuation  for  given  frequencies44 


One  can  see  from  Figure  1 1  that  the  operating  frequency  of  300  GHz,  while  on  the  rise  of  a  relative  lull  of 
atmospheric  attenuation  still  misses  the  peak  of  attenuation.  Thus,  this  frequency  lends  well  to  be  used  for 
long-range  transmission  in  order  to  interrogate  a  solid  rocket  exhaust  plume. 

Another  consideration  lies  in  the  operation  of  the  antenna.  In  general  there  are  two  modes  to 
operate  in:  the  far-field  and  the  near-field.  Compared  with  the  far-field  mode  of  operation,  the  near-field  is 
a  complicated  area  to  understand.  While  there  has  been  work  done  to  characterize  the  near-field,  problems 
that  can  be  encountered  are  the  stored  energy  of  a  signal  and  reactive  power.45  Thus,  to  avoid  the 
complications  inherent  to  operating  in  the  near-field,  this  study  conducted  all  of  its  theoretical  and 
experimental  work  in  the  far-field.  To  operate  in  the  far- field,  certain  criteria  must  be  met  which  puts 
restrictions  on  minimum  distances  of  operation  based  on  operating  wavelength  and  antenna  dimensions. 
The  criteria  for  far-field  that  were  used  for  this  study  are  defined  by  Stutzman46  as 


R  »  A 


R  »  Ddish 


R  > 


2D dish 
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One  can  see  from  these  equations  that  depending  on  the  setup  a  different  criterion  will  dominate.  Later  in 
this  study,  it  will  be  shown  that  the  proof-of-concept  experiment  is  dominated  by  the  second  condition, 
while  the  notional  terahertz  diagnostic  instrument  must  satisfy  the  third  condition. 

2.  Parabolic  Reflector  Antennas 

Reflector  antennas  have  been  studied  for  many  years.  They  are  generally  derived  from  the  conic 
sections:  circle,  ellipse,  parabola,  hyperbola,  and  straight  line.  The  first  reflector  antenna  was  a  cylindrical 
parabolic  reflector  created  by  Heinrich  Hertz  in  1888.  Antennas  were  first  used  for  experiments,  but  then 
as  use  of  technology  in  warfare  increased,  they  were  developed  for  radar  and  communications.  Since  the 
dawn  of  the  space  age,  reflector  antennas  have  made  large  jumps  in  terms  of  size,  gain,  and  application  in 
order  to  keep  up  with  the  demand  to  track  and  communicate  with  space  vehicles.  Today,  some  of  the  most 
prominent  reflector  antennas  include  the  Arecibo  1000-ft  spherical  reflector,  the  ATS-6  unfurlable  30-ft 
reflector,  and  the  famed  Effelsberg  100-meter  radio  telescope  which  has  a  maximum  realizable  gain  of  83 
dB.47  These  antennas  are  pictured  below  in  Figure  12.  Seeing  the  capabilities  of  such  reflector  antennas,  it 
was  decided  to  use  a  parabolic  reflector  antenna  in  the  design  of  the  notional  terahertz  diagnostic 
instrument. 
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(c) 

Figure  12:  Several  large  reflector  antennas  (a)  Effelsberg  100-meter  radio  telescope  (b)  ATS-6 
unfurlable  300-ft  reflector  (c)  Arecibo  1000-ft  spherical  reflector47 


In  general,  parabolic  reflector  antennas  operate  by  using  a  smaller  low  gain  antenna  as  a  feed 
antenna  which  generates  the  signal  inside  the  reflector.  The  emitted  signal  is  then  reflected  off  of  the 
antenna  dish  and  propagated  outward  from  the  antenna.  The  factor  of  most  concern  for  application  in  the 
notional  instrument  is  the  antenna  gain.  Free  space  loss,  as  discussed  previously,  poses  the  largest  threat  to 
the  operation  of  this  system.  If  too  much  attenuation  occurs  on  the  way  to  and  from  the  solid  rocket 
exhaust  plume,  the  scattering  information  will  be  lost  in  atmospheric  noise  and  be  useless  in  terns  of 
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performing  a  diagnostic  of  the  particles  in  the  plume.  The  gain  for  a  parabolic  reflector  can  be  calculated 
using  the  following  equation.48 
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An  (nDjish\( 
A2(  4  j( 


a2y3 
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The  information  inside  the  brackets  is  generally  held  as  the  gain  factor  of  the  antenna.  For  the  purposes  of 
designing  the  notional  terahertz  diagnostic  instrument,  a  simplified  version  of  this  equation  was  used.  It 
depends  only  on  the  diameter  of  the  dish  and  the  frequency  of  the  transmitted  light.  It  assumes  a  gain 
factor  of  0.55  and  is  used  in  the  development  of  the  notional  instrument  for  its  ease  of  use.49 

G  =  17.8  +  20  log Ddish  +  201og  /  (20) 

This  equation  combined  with  the  restrictions  for  far-field  operation  set  the  basis  for  the  design  of  the 
notional  terahertz  solid  rocket  exhaust  plume  diagnostic  tool. 

3.  Horn  Antennas 

The  proof-of-concept  experiment  was  conducted  in  an  indoor  chamber  that  would  not  allow  for 
the  use  of  a  parabolic  antenna  while  satisfying  the  far-field  criteria.  Therefore,  a  different  antenna  type  was 
selected:  the  horn  antenna,  which  is  a  variety  of  aperture  antenna.  There  are  three  types  of  rectangular  horn 
antennas:  (1)  E-plane  sectoral,  which  is  flared  only  in  the  E-plane,  (2)  H-plane  sectoral,  which  is  flared 
only  in  the  H-plane,  and  (3)  a  pyramidal  horn,  which  is  flared  in  both  the  E  and  H  planes.  The  difference  in 
the  horn  varieties  can  be  seen  in  Figure  13. 


Figure  13:  Horn  antenna  types:  (a)  H-plane  sectoral  (b)  E-plane  sectoral  (c)  Pyramidal50 
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The  proof  of  concept  experiment  utilizes  only  the  pyramidal  horn  antenna,  thus  this  section  will  only  look 
at  the  theory  for  that  specific  type. 

In  general,  a  horn  antenna  provides  a  smooth  transition  for  electromagnetic  radiation  from  a 
waveguide  into  free  space  or  vice  versa.  It  operates  in  much  the  same  manner  as  a  megaphone  amplifies 
the  sound  of  someone’s  voice.  Again,  the  important  issue  for  this  study  is  antenna  gain.  However,  for  the 
proof-of-concept  experiment  the  size  of  the  horn  must  also  be  taken  into  account  so  that  the  far-field 
criteria  are  not  violated.  Horn  antennas  are  characterized  by  their  H-plane  and  E-plane  patterns  which  are 
derived  from  the  antenna’s  gain.  The  gain  can  be  calculated  using  the  following  equation  assuming  an 
aperture  efficiency  of  50%. 50 


G  = 


1  47T 
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In  this  equation  A  and  B  represent  the  horizontal  and  vertical  dimensions  of  the  aperture  respectively. 
These  two  parameters  are  also  used  to  determine  the  far-field  restrictions  that  are  placed  on  the  experiment. 
The  largest  dimension  of  the  horn  antenna  is  the  diagonal  of  the  aperture  which  can  be  easily  found  using 
the  horizontal  and  vertical  dimensions.  This  dimension  is  then  applied  to  the  far- field  criteria  in  Eq.  (18)  to 
determine  the  minimum  distance  for  both  the  transmitter  and  the  receiver. 

C.  Rayleigh  Approach 

The  intensity  of  light  scattered  from  a  particle  is  a  function  of  the  incident  intensity  of  light, 
particle  radius,  distance,  and  wavelength  as  shown  in  Eq.  (15).  Eq.  (15)  represents  the  total  intensity  of 
light  scattered.  In  general,  it  is  not  feasible  to  collect  all  of  the  scattered  light.  In  fact,  only  a  small  portion 
of  this  light  can  be  measured  through  the  collection  aperture.  In  order  to  better  analyze  the  situation,  one 
can  define  a  differential  cross-section. 
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From  Miles,  Lempert,  and  Forkey40,  the  total  collected  power  becomes  the  integral  of  the  number  of 
particles  multiplied  by  the  scattering  differential  cross-section,  the  incident  intensity,  and  the  scattering 
efficiency  over  the  collection  solid  angle.  Namely, 
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The  differential  solid  angle  can  be  represented  as 
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which  when  applied  to  the  integral  for  detected  power  becomes 
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Evaluating  the  integral  in  Eq.  (25)  for  the  power  yields  the  equation 

97r2E2  (m2  —  1\2  ( ( 91  —  02)(cos(301)  —  cos(302)  —  cos (90x)  +  cos  (9 02))\  (26) 
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For  a  comparison  between  Mie  theory  and  Rayleigh  theory,  Cox,  DeWeerd,  and  Linden51 
performed  an  experiment  to  measure  the  Mie  and  Rayleigh  total  scattering  cross-section  for  several 
particles.  The  experiment  used  three  different  particle  radii  and  a  frequency  sweeping  source.  They 
confirmed  that  while  the  criterion  for  Rayleigh  scattering  is  met,  the  two  theories  produce  the  same  results. 

In  a  particular  experimental  setup,  certain  values  can  be  controlled  and  are  known.  The  values  for 
incident  intensity,  wavelength,  refractive  index,  and  collection  angle  are  generally  known.  The  experiment 
will  typically  yield  a  value  for  the  scattered  power.  Unfortunately,  this  leaves  two  unknown  values  and 
only  one  equation.  Those  values  are  particle  number  density  and  particle  radius.  In  past  experiments 
utilizing  Rayleigh  scattering,  the  number  density  was  either  known  or  reasonably  approximated.  For 
example,  Zhang  et  al 35  measured  the  size  and  number  density  of  glycogen  molecules.  However,  in  their 
experiment,  the  size  of  glycogen  was  relatively  well  known  as  well  as  the  concentrations  of  the  solutions 
they  used  which  gives  a  relation  to  number  density.  Also,  An-Le  et  al30  measured  the  size  of  argon  clusters 
using  Rayleigh  scattering.  However,  they  were  able  to  relate  both  size  and  number  to  their  reservoir 
pressure.  Unfortunately,  these  values  are  not  generally  known  for  alumina  particles  in  solid  rocket  plumes. 
There  have  been  studies  conducted  to  determine  particle  size  based  on  the  propellant  characteristics.  O.  B. 
Kovalev52  discussed  the  formation  of  aluminum  oxide  in  rocket  plumes  and  generated  an  analytical 
prediction  method  for  the  size  of  the  particles.  However,  the  particle  number  density  and  size  in  a  typical 
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rocket  plume  is  generally  unknown.  Therefore,  this  technique  may  be  restricted  to  measuring  the  product 
of  NV2  unless  other  methods  can  be  developed  in  the  future.  One  method  would  involve  sweeping 
frequencies.  This  technique  would  start  at  the  smallest  frequency  in  the  terahertz  regime  and  increase  by 
discrete  amounts.  This  increase  would  remove  the  scattering  signal  from  the  larger  particles  because  they 
would  no  longer  satisfy  the  Rayleigh  criterion.  These  particles  would  be  in  the  Mie  regime  where  the  off 
incident  direction  scattering  is  negligible.  The  change  in  scattered  signal  return  can  be  measured  and 
analyzed  to  calculate  a  number  density  for  a  certain  group  of  particle  radii.  The  process  could  be  repeated 
to  account  for  all  particle  sizes  of  interest.  Also,  using  interferometry  could  lead  to  measuring  the  particle 
velocity.53 


In  order  to  calculate  the  actual  power  received  from  a  scattering  center,  the  transmission  losses  and 
gains  from  antennas  must  be  taken  into  account.  Friis  transmission  formula  is  used  to  determine  the  power 
received  based  on  the  space  loss  from  Eq.  (17)  and  the  gains  of  the  antenna. 
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Using  this  formula  with  Eq.  (25),  the  power  ratio  for  an  actual  setup  can  be  written  as 
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CHAPTER  6 


EXPERIMENTAL  SETUP  AND  PREDICTED  RESULTS 
A.  Experimental  Setup 

The  proof-of-concept  experiment  was  designed  to  be  a  scaled  version  of  the  notional  terahertz 
diagnostic  instalment  working  with  typical  solid  particles  in  a  solid  rocket  plume.  Therefore,  the  sphere 
sizes  were  scaled  with  the  experimental  wavelength  (12.24  cm)  in  order  to  meet  the  Rayleigh  scattering 
criterion,  reproducing  the  plume  environment  as  far  as  particle  size  to  wavelength  ratio  is  concerned.  The 
Electrical  and  Computer  Engineering  Department  at  the  University  of  Colorado  at  Colorado  Springs  is 
currently  developing  the  terahertz  source  that  will  eventually  be  used  in  this  project.  This  source  is 
designed  to  operate  at  300  GHz  with  an  output  power  of  1  kW  with  minimal  frequency  sweeping 
capabilities.  Assuming  that  the  notional  terahertz  diagnostic  instillment  would  have  a  non-frequency 
sweeping  source,  the  wavelength  of  1  mm  would  dictate  a  maximum  particle  radius  of  100  microns  in  order 
to  stay  within  Rayleigh  theory  constraints.  The  proof-of-concept  experiment  utilized  microwaves  at  a 
frequency  of  2.45  GHz  and  wavelength  of  12.24  cm.  This  source  had  a  fixed  frequency  which  set  the 
maximum  particle  radii  at  12.24  mm.  The  experiment  looked  at  particles  whose  radius  ranged  from  4.76 
mm  to  9.53  mm. 

All  of  the  experiments  were  conducted  inside  a  microwave  anechoic  chamber  at  the  University  of 
Colorado  at  Colorado  Springs.  The  setup  consisted  of  a  power  source  transmitted  through  a  horn  antenna. 
The  scattered  energy  from  the  particle  was  collected  in  a  second  horn  antenna.  The  power  source  used  in 
the  experiment  is  a  National  Electronics  2  kW  magnetron  source  which  was  coupled  with  a  National 
Electronics  Isolator,  a  waveguide  measurement  piece,  and  an  Agilent  Tuner.  The  power  source  along  with 
attached  horn  antenna  can  be  seen  in  Figure  14. 


Distribution  A:  Approved  for  public  release.  Distribution  unlimited. 


35 


Figure  14:  2  kW  magnetron  source  and  horn  antenna 


The  forward  power  from  the  2  kW  magnetron  source  was  measured  using  a  HP  power  sensor 
connected  to  an  Agilent  power  meter.  The  measured  data  was  recorded  using  Lab  View  8.2.1. 
Microwaves  were  passed  through  coaxial  cables  and  emitted  into  the  anechoic  chamber  from  a  10  dB 
standard  gain  horn.  The  microwaves  traveled  approximately  2.44  m  to  the  sphere  as  shown  in  Figure  15 
satisfying  the  far-field  requirement  defined  by  Eq.  (18). 


Figure  15:  Distance  traveled  by  microwaves  from  transmitter 
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The  spheres  were  made  out  of  three  different  materials:  2017-T4  aluminum,  brass,  and  stainless 
steel  302.  The  spheres  had  several  different  radii  which  are  9.53  mm,  7.94  mm,  6.35  mm,  and  4.76  mm. 
They  were  supported  with  6.35  mm  outer  diameter  (4  mm  inner  diameter)  quartz  tubes  that  were 
approximately  14  cm  apart  for  the  double  sphere  setup.  The  4.76  mm  radius  aluminum  spheres  can  be  seen 
on  the  quartz  support  tubes  in  Figure  16  in  both  a  single  and  a  double  sphere  setup. 


Scattering  Spheres  Anechoic  Cones  Quartz  Rods 


(A)  CB) 

Figure  16:  Sphere  configuration  (A)  Single  and  (B)  Double 


The  scattered  microwaves  were  collected  in  a  14.7  dB  standard  gain  horn.  The  received  power 
was  passed  through  coaxial  cables  out  the  anechoic  chamber  to  an  HP  power  sensor  which  was  connected 
to  the  Agilent  power  meter.  The  receiving  horn  can  be  seen  in  Figure  17  connected  to  the  power  sensor. 
Figure  18  and  Figure  19  show  the  entire  setup  from  the  perspective  of  the  spheres  and  an  observer.  Figure 
20  provides  a  general  overhead  view  of  the  experimental  setup. 

Power  Sensor  Receiver  Horn 


Figure  17:  Receiving  horn  for  scattered  microwaves 
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Figure  18:  Experimental  setup  from  observer  perspective 


Transmitting  Horn 


Receiver  Horn 


Figure  19:  Experimental  setup  from  sphere  perspective 


Figure  20:  Overhead  view  of  experimental  setup 
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The  concept  of  using  microwave  horns  as  both  the  transmitting  and  receiving  antennas  is 
consistent  with  Zhang,  Shneider,  and  Miles'8  who  used  a  similar  setup  to  measure  Rayleigh  scattering  using 
microwaves  from  resonance-enhanced  multiphoton  ionization  in  argon.  They  used  a  10  mW  source  at  a 
frequency  of  12.6  GHz  and  collected  the  scattered  signal  in  a  microwave  homodyne  receiving  system. 
Shneider  and  Miles36  also  performed  microwave  Rayleigh  scattering.  Using  frequencies  of  3  and  6  GHz, 
they  performed  plasma  diagnostics  by  treating  the  plasma  as  a  point  dipole  source  and  collecting  the 
Rayleigh  scattered  light  from  the  plasma. 

The  experiments  performed  in  this  study  were  repeated  for  each  of  the  sphere  sizes  and  materials. 
The  input  power  from  the  microwave  source  ranged  from  approximately  200  Watts  to  2000  Watts.  The 
experiment  was  also  run  with  only  the  quartz  rod  in  the  chamber.  The  collected  power  from  the  “no 
sphere”  run  was  subtracted  from  each  of  the  recorded  tests  with  a  sphere  in  the  chamber  as  a  zero  signal 
reference.  The  difference  was  taken  as  the  scattered  power  and  was  used  in  comparison  with  the  other 
tests.  Eight  data  points  were  taken  which  represent  the  average  of  multiple  individual  tests.  The  eight 
points  corresponded  to  input  powers  that  are  approximately  200  Watts  to  2000  Watts  with  a  200  Watt 
interval. 

B.  Predicted  Experimental  Results 

The  aforementioned  proof-of-concept  experimental  setup  parameters  were  input  into  a  MATLAB 
code  in  order  to  generate  predicted  results  for  the  tests  conducted.  The  known  parameters  for  the  proof-of- 
concept  experiment  include  the  gains  of  the  horn  antennas,  the  geometry  of  the  experiment,  and  the  free 
space  loss  for  transmission.  The  linear  gain  for  the  transmitting  and  receiving  antennas  were  measured  to 
be  10  and  29.51  respectively.  The  free  space  loss  for  the  transmission  of  the  input  power  is  1.59xl0_s. 
The  geometry  determines  the  bounds  of  the  integral  term  from  Eq.  (28)  which  yields  a  value  of  0.021. 
Inputting  these  parameters  into  Eq.  (28),  it  can  be  simplified  to 


-f  =  18. 798917 IV  vr4 


m2  —  1 
m2  +  2 


(29) 


Eq.  (29)  provides  the  basis  for  the  Matlab  code  which  can  be  seen  in  Appendix  B. 

The  two  unknown  material  properties  in  Eq.  (29)  are  the  index  of  refraction  and  the  scattering 
efficiency.  Reference  54  investigated  values  for  the  index  of  refraction  based  on  different  wavelengths  of 
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light  for  aluminum.  Figure  21  shows  the  results  for  the  relation  between  wavelength  and  index  of 
refraction.  It  is  obvious  that  as  wavelength  continues  to  increase  the  index  of  refraction  also  continues  to 
increase. 


Figure  21:  Relation  between  refractive  index  and  wavelength  for  aluminum 

The  relation  in  Figure  21  agrees  with  Reference  55  which  states  that  as  the  frequency  decreases, 
the  index  of  refraction  for  metals  increases  to  infinity.  Figure  22  and  Figure  23  show  the  expected  trends 
for  the  proof-of-concept  experiment.  In  general,  the  index  of  refraction  at  microwave  and  terahertz 
frequencies  for  various  materials  is  unknown.  A  thorough  review  of  the  existing  literature  yielded  no 
conclusive  results.  Therefore,  an  assumed  index  of  refraction  of  2  and  a  scattering  efficiency  of  100%  were 
used  in  Figure  22  and  Figure  23  so  that  general  trends  can  be  seen.  Figure  22  shows  the  linear  trend  from 
Eq.  (29)  between  the  source  power  and  the  scattered  power  as  well  as  the  increase  in  scattered  power 
associated  with  an  increase  in  the  number  of  scattering  centers  for  a  sphere  with  a  radius  of  9.53  mm. 
Figure  23  depicts  the  fourth  order  relationship  between  the  radius  of  the  scattering  sphere  and  the  scattered 
power  as  defined  by  Eq.  (29).  Figure  24  highlights  the  role  that  the  index  of  refraction  plays  in  the 
scattered  power  collected  from  a  sphere  with  a  radius  of  9.53  mm  and  a  scattering  efficiency  of  100%.  It  is 
clear  that  the  scattered  power  is  not  substantially  affected  by  the  index  of  refraction. 
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Figure  22:  Scattered  power  versus  source  power  for  generic  sphere  with  r  =  9.53  mm  and  m  =  2 


Figure  23:  Scattered  power  versus  sphere  radius  with  m  =  2 


Figure  24:  Scattered  power  versus  source  power  for  a  9.53  mm  radius  sphere  as  a  function  of  the 

refractive  index 
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The  effect  of  the  scattering  efficiency  must  also  be  taken  into  account  for  the  experiment  since  it  is 
also  not  generally  known  at  the  wavelengths  of  interest  in  this  study.  While  Figure  24  shows  that  an 
increase  in  the  refractive  index  leads  to  an  increase  in  the  scattered  power,  this  increase  can  be  offset  by  a 
decrease  in  the  scattering  efficiency  of  a  particle.  Figure  25  shows  this  phenomenon  by  holding  the  index 
of  refraction  fixed  at  109  and  varying  the  scattering  efficiency. 
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Figure  25:  Scattered  power  versus  source  power  for  a  fixed  refractive  index  of  m  =  109 
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CHAPTER  7 


RESULTS  AND  DISCUSSION  OF  PROOF  OF  CONCEPT  EXPERIMENT 
A.  Single  Sphere  Results 

The  first  tests  done  for  the  proof-of-concept  used  only  a  single  sphere  as  the  scattering  center. 
Only  two  materials  were  used  for  these  tests:  aluminum  and  brass.  The  results  for  these  tests  can  be  seen  in 
Figure  26  and  Figure  27. 


r  =  9.53  mm 
r  =  4.76  mm 


Figure  26:  Scattered  power  versus  source  power  for  one  aluminum  sphere 


r  =  9.53  mm 
r  =  7.94  mm 
r  =  6.35  mm 
r  =  4.76  mm 


Figure  27:  Scattered  power  versus  source  power  for  one  brass  sphere 

In  Figure  26  and  Figure  27,  the  linear  trend,  as  defined  by  Eq.  (29),  is  apparent.  Flowever,  the 
error  bars  for  the  measurements  are  very  large.  The  error  stems  from  the  small  signal  to  noise  ratio 
encountered  in  these  tests.  In  Figure  26,  the  average  signal  to  noise  ratio  is  1.06  for  the  9.53  mm  radius 
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sphere  and  1.004  for  the  4.76  mm  radius  sphere.  In  Figure  27,  the  signal  to  noise  ratio  ranges  from  1.005  to 
1.031.  For  accurate  measurements,  a  larger  signal  to  noise  ratio  is  needed.  In  order  to  improve  the  signal 
to  noise  ratio,  the  remainder  of  the  tests  used  two  spheres  which  acts  to  approximately  double  the  scattered 
signal  according  to  Eq.  (29). 

B.  Double  Sphere  Results 

The  next  set  of  tests  were  conducted  with  two  spheres  for  reasons  mentioned  above.  Each  test 
data  point  was  repeated  five  times  and  the  average  of  these  tests  is  reported.  The  error  bars  on  the  graphs 
represent  the  standard  deviation  for  the  group  of  tests.  Each  different  material  and  sphere  size  was  used  in 
these  tests.  The  results  for  these  runs  can  be  seen  in  Figure  28,  Figure  29,  and  Figure  30. 
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Figure  28:  Scattered  power  versus  source  power  for  two  aluminum  spheres 


The  predicted  line  in  Figure  28  was  generated  using  Eq.  (29)  and  assumed  values  for  the  unknown 
index  of  refraction  and  scattering  efficiency.  The  scattering  efficiency  was  assumed  to  be  100%,  and  an 
index  of  refraction  of  3.4  provided  the  best  fit  to  the  experimental  results  for  the  9.53  mm  radius  sphere. 
Similarly,  values  for  the  index  of  refraction  were  chosen  to  produce  best  fit  prediction  lines  for  Figure  29 
and  Figure  30.  The  value  for  brass  is  6.0,  and  the  value  for  stainless  steel  is  1.85.  While  the  true  values 
remain  unknown,  the  determination  of  these  parameters  is  beyond  the  scope  of  this  study.  However,  for  a 
terahertz  diagnostic  instrument,  the  actual  values  for  alumina  will  need  to  be  determined  for  the  operating 
wavelengths  and  temperatures. 


Distribution  A:  Approved  for  public  release.  Distribution  unlimited. 


44 


=L 


Oi 

g 

o 

Q. 

T5 

<U 


(t! 

u 

l/l 


500 


1000  1500 

Source  Power  (W) 


2000 


2500 


Figure  29:  Scattered  power  versus  source  power  for  two  brass  spheres 


Figure  30:  Scattered  power  versus  source  power  for  two  stainless  steel  spheres 


In  each  of  the  tests  there  is  a  linear  trend  for  the  experiment,  which  is  expected  based  on  Eq.  (29). 
Comparing  these  figures  to  that  of  the  single  sphere  data,  one  can  clearly  see  the  increase  in  the  scattered 
signal  along  with  the  decrease  in  the  error  bars.  The  signal  to  ratio  for  these  tests,  however,  is  still  low 
ranging  from  1.009  to  1.197.  While  it  may  be  the  case  for  each  of  the  materials  in  the  experiment  that  the 
index  of  refraction  is  much  higher  than  the  assumed  value,  the  scattering  efficiency  can  offset  the  increase 
in  scattered  power.  Using  109  as  the  refractive  index  for  each  material,  a  scattering  efficiency  of  60.65% 
for  aluminum  would  create  the  same  predicted  trend  shown  in  the  results  based  on  Eq.  (29).  Likewise,  a 
scattering  efficiency  of  84.83%  readjusts  the  brass  predicted  trend,  and  an  efficiency  of  19.96%  for 
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stainless  steel  satisfies  the  trend  difference.  The  reality  for  the  materials  is  probably  a  compromise  between 
these  two  values  in  order  to  generate  the  results  from  the  proof-of-concept  experiment.  It  is  obvious  that 
for  any  experiment  or  diagnostic,  accurate  values  for  these  parameters  are  necessary. 

C.  Sphere  Number  Comparison 

Figure  31  shows  the  ratio  of  the  scattered  power  between  two  spheres  and  a  single  sphere  for  the 
9.53  mm  radius  sphere  of  aluminum. 
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Figure  31:  Aluminum  sphere  number  comparison 

In  Figure  31,  the  highest  returned  scattered  power  for  the  double  sphere  set  up  is  372.6  pW,  and  the 
single  sphere  setup  maximum  is  178.3  pW.  Within  the  errors  bars  of  the  two  different  setups,  there  is  a 
factor  of  two  increase  in  the  scattered  power  as  expected  by  Eq.  (29);  therefore,  doubling  the  scattering 
centers  leads  to  an  effective  doubling  of  the  collected  Rayleigh  scattered  signal.  The  low  signal  to  noise 
ratio  is  the  source  of  the  fluctuations  above  the  ratio  of  2. 

D.  Sphere  Size  Comparison 

The  next  comparison  for  the  collected  data  was  between  the  different  sphere  sizes.  According  to 
Eq.  (29)  the  relation  between  the  different  sphere  radii  should  be  fourth  order.  Thus  for  Figure  32,  Figure 
33,  and  Figure  34  the  different  sphere  sizes  were  compared  within  the  same  material.  For  each  of  the 
graphs,  the  highest  output  power  from  the  source  was  used  in  order  to  get  the  highest  signal  to  noise  ratio. 
Also,  only  the  double  sphere  data  was  used  in  creating  the  following  figures. 
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Figure  32:  Aluminum  sphere  size  comparison 
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Figure  33:  Brass  sphere  size  comparison 


Figure  34:  Stainless  Steel  sphere  size  comparison 
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In  Figure  32,  Figure  33,  and  Figure  34,  the  fourth  order  relation  between  radius  and  scattered  power 
agrees  with  Eq.  (29).  The  predicted  line  in  these  figures  was  generated  using  the  best  fit  index  of  refraction 
found  earlier  in  this  chapter.  Those  values  are  3.4,  6.0,  and  1.85  for  aluminum,  brass,  and  stainless  steel 
respectively.  The  largest  sphere  radius  was  used  to  determine  the  best  fit  index  of  refraction  because  it 
yielded  the  highest  signal  to  noise  ratio  and,  therefore,  accuracy.  However,  the  fit  for  the  largest  sphere 
radius  is  not  perfect,  and  was  modified  in  order  to  improve  the  fit  for  the  smaller  radii.  It  can  also  be  seen 
that  even  with  a  modified  index  of  refraction  the  theory  and  experiment  are  not  a  perfect  match  especially 
at  smaller  sphere  radii.  The  discrepancy  may  be  due  to  uncertainties  in  the  index  of  refraction  as  well  as  in 
the  scattering  efficiency.  Also,  the  small  signal  to  noise  ratio  is  still  a  factor  in  the  experimental  error. 

E.  Sphere  Material  Comparison 

The  final  comparison  for  the  data  from  the  proof-of-concept  scattering  experiment  was  between 
the  three  different  sphere  materials.  The  difference  between  the  materials  should  be  based  solely  on  the 
difference  in  the  index  of  refraction  and  scattering  efficiencies  between  the  different  materials.  Figure  35 
and  Figure  36  show  the  scattered  power  versus  source  power  curves  for  the  sphere  sizes  of  9.53  mm  and 
4.76  mm. 
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Figure  35:  Sphere  material  comparison  for  9.53  mm  radius  sphere 
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Figure  36:  Sphere  material  comparison  for  4.76  mm  radius  sphere 


It  can  be  seen  in  Figure  35  and  Figure  36  that  the  index  of  refraction  and  the  scattering  efficiency  directly 
affect  the  returned  scattered  power  as  Eq.  (29)  indicates. 

F.  Experiment  Discussion 

The  results  of  the  proof-of-concept  experiment  show  that  Rayleigh  scattering  can  be  used  to 
determine  particle  size  when  the  number  density  of  the  particles  is  known.  The  data  is  in  good  agreement 
with  the  trends  set  forth  by  the  theoretical  calculations  for  the  experiment.  However,  there  are  some 
differences  between  theory  and  experiment  and  sources  of  error  that  need  to  be  addressed. 

1.  Sources  of  Error 

There  were  several  sources  of  error  that  arose  during  the  process  of  conducting  the  proof-of- 
concept  experiment.  One  of  the  most  prominent  sources  of  error  between  the  theory  and  the  experiment 
was  the  small  signal  to  noise  ratios  generated  during  the  tests.  The  maximum  ratio  from  the  tests  was  1.18 
which  is  well  below  the  threshold  for  accurate  data.  The  low  ratio  led  to  large  error  bars  in  the 
experiments.  The  large  amount  of  noise  was  mostly  generated  by  reflection  inside  the  anechoic  chamber. 
The  reflection  led  to  milli-Watts  of  measured  power  for  the  zero  traces  of  the  experiment.  Another  large 
source  of  error  was  the  uncertainty  in  the  index  of  refraction  at  microwave  frequencies  and  the  scattering 
efficiency  of  the  materials  used  in  the  experiment.  These  parameters  have  a  strong  influence  on  the 
predicted  amount  of  scattered  power  for  each  of  the  tests. 
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Other  sources  of  error  include  the  power  sensors  used  to  measure  both  the  transmitted  and 
scattered  signals.  These  sensors  have  a  quoted  measurement  noise  of  <110nW  which  is  0.16%  of  the 
lowest  measured  signal.  Another  factor  would  be  minor  fluctuations  in  the  transmitted  signal.  For  the 
purposes  of  ease  of  measurement,  the  magnetron  source  was  allowed  to  fluctuate  by  +/-0.64  W.  Any 
movement  of  the  anechoic  cones  in  the  chamber  could  also  have  an  effect  on  the  consistency  between 
different  test  runs.  Another  experimental  setup  factor  that  would  play  into  the  error  bars  is  the  quartz  rods. 
While  mostly  microwave  transparent,  they  still  contributed  slightly  to  the  measured  powers.  Again  if  the 
rods’  position  is  moved  slightly  the  reflected  signal  could  change  leading  to  inconsistencies  between  the 
individual  runs.  Finally,  any  discrepancies  on  the  spheres  themselves  would  lead  to  a  change  in  signal 
based  on  their  orientation  on  the  quartz  rods.  The  spheres  were  not  perfectly  spherical,  and  they  were  not 
set  on  the  rods  the  same  every  time. 

2.  Differences  Between  Theory >  and  Results 

The  main  difference  between  theory  and  results  was  in  the  magnitude  of  the  returned  scattered 
power.  The  experimental  values  for  scattered  power  were  only  slightly  different  than  the  predicted  results 
for  the  larger  sphere  radii.  Flowever,  the  percent  differences  grew  for  the  smaller  sphere  radii.  Table  5 
shows  the  percent  difference  between  the  predicted  results  and  the  experimental  data  for  the  double  sphere 
tests. 


Refractive  Index 

%  Difference 

3.40 

9.53  mm  Al 

-5.45 

4.76  mm  Al 

113.12 

6.00 

9.53  mm  Br 

-4.40 

7.94  mm  Br 

15.06 

6.35  mm  Br 

49.44 

4.76  mm  Br 

88.11 

1.85 

9.53  mm  St 

4.18 

7.94  mm  St 

21.50 

6.35  mm  St 

44.13 

4.76  mm  St 

159.06 

Table  5:  Differences  between  predicted  and  experimental  results 
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It  is  important  to  note  again  that  the  scattering  efficiency  for  the  indices  of  refraction  was  assumed 
to  be  100%.  It  is  most  likely  a  tradeoff  between  the  actual  index  of  refraction  with  a  lower  scattering 
efficiency.  The  increase  in  percent  difference  as  the  sphere  size  decreases  is  most  likely  a  result  of  the 
uncertainty  in  this  trade  off  and  the  low  signal  to  noise  ratio  that  occurs  at  the  smaller  radii.  Despite  these 
differences,  the  experimental  results  still  show  the  capability  of  calculating  the  particle  radius  based  on  the 
scattered  data.  It  is  clear  that  each  of  the  sphere  materials  and  sizes  had  a  unique  scattered  power  which 
would  enable  identifying  a  particle  based  on  its  scattered  signal. 
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CHAPTER  8 


DEVELOPMENT  OF  NOTIONAL  TERAHERTZ  INSTRUMENT 


A.  Notional  Terahertz  Instrument 

The  proof-of-concept  experiment  demonstrated  the  capability  of  measuring  particle  size  based  on 
the  amount  of  scattered  light  received  if  the  number  density  of  the  particles  is  known.  In  order  to  extend 
the  principles  explored  in  the  experiment  to  the  diagnostic  of  a  solid  rocket  exhaust  plume,  the  Rayleigh 
theory  and  far-field  criteria  must  be  revisited. 

For  an  actual  experimental  setup  involving  terahertz  radiation  and  Rayleigh  scattering,  one  must 
consider  the  challenges  of  transmitting  across  large  distances  while  scattering  off  particles  inside  the  rocket 
plume.  For  the  predictions  of  the  required  power  from  a  source,  certain  assumptions  were  made.  For 
Rayleigh  theory  to  apply,  the  wavelength  was  ten  times  greater  than  the  particle  radius.  Also,  the  criteria 
for  far- field  transmission  were  used  as  defined  in  Eq.  (18)  by  Stutzman.46  After  satisfying  these  criteria,  the 
Friis  transmission  formula42  from  Eq.  (27)  was  used  to  calculate  the  ratio  between  the  power  at  the 
scattering  particle  and  the  power  transmitted.  The  notional  setup  uses  parabolic  antennas  to  transmit  and 
receive.  The  gains  for  the  antennas  are  based  on  frequency  and  the  radius  of  the  dish  and  can  be  calculated 
with  Eq.  (20). 49  With  this  setup,  the  integral  bounds  of  the  scattering  cross-section  can  be  defined  in  terms 
of  the  distance  to  the  particle  and  the  radius  of  the  antenna  dish;  namely, 


The  final  ratio  between  received  power  and  transmitted  power  can  be  written  as 
P,  /  A  \2  9n2V2  ( n 2  -  lV  f  f 

W  -(rfTz)  J  J  sin 


-  GtGr 


(32) 


9±  (p1 


Figure  37  shows  the  effects  of  distance  and  antenna  size  on  the  power  ratio  Eq.  (32)  for  a  single 
100  micron  scattering  particle.  The  graphs  show  that  the  while  a  large  antenna  provides  more  gain,  the 
effects  of  increased  distance,  due  to  the  far-field  restrictions,  are  a  greater  detriment  to  the  power  ratio. 
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Overall,  the  relation  is  inversely  proportional  to  the  fourth  power  of  the  range.  Figure  38  shows  the  relation 
between  antenna  radius  and  the  power  ratio  assuming  that  the  range  is  at  the  minimum  distance  while  still 
satisfying  the  far-field  criteria.  Again  it  shows  the  inverse  proportionality  to  the  fourth  power  of  the 
distance  to  the  particle.  The  Matlab  code  used  to  generate  these  plots  can  be  seen  in  Appendix  C. 


12.00 


Antenna 

. 0.75 

- 1  m 

—  .  .1,25  m 
-  .  1.5  m 


Range  to  Particle  (m) 


Figure  37:  Relation  between  range  to  particle  and  power  ratio  with  multiple  antenna  radii 


Figure  38:  Relation  between  antenna  radius  and  power  ratio  assuming  minimum  range  to  particle 

From  Figure  37  and  Figure  38,  one  can  see  that,  for  a  notional  terahertz  diagnostic  instrument,  a 
high  power  source  is  required  with  high  gain  antennas.  These  figures  show  that  the  biggest  detriment  to 
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scattered  power  is  distance.  Also,  the  power  ratio  for  a  single  particle  is  very  small.  However,  the  number 
density  for  solid  rocket  plumes  can  be  calculated  for  a  Castor  120  solid  rocket  motor  to  be  on  the  order  of 
10s  to  1010  particles  per  cubic  meter  using  a  mass  flow  of  3.54  kg/s  and  a  particle  mass  percentage  of 
32%. 56  This  number  density  woidd  allow  a  1  kW  source  using  1.5  m  radius  antennas  to  measure  scattered 
powers  in  the  pW  range  at  a  distance  of  20  km.  Figure  39  depicts  what  the  notional  setup  may  look  like. 


Figure  39:  Notional  terahertz  diagnostic  instrument 


B.  Complications  of  Implementation  for  Notional  Terahertz  Instrument 

The  two  largest  complications  arising  for  the  implementation  of  the  notional  terahertz  diagnostic 
instrument  are  the  power  source  and  the  antenna.  As  discussed  in  Chapter  2,  there  is  a  current  lack  of  high 
power  sources  that  operate  in  the  terahertz  regime.  While  the  1  kW  source  that  was  also  previously 
discussed  would  prove  beneficial  to  the  terahertz  diagnostic  tool,  it  still  would  only  be  able  to  generate 
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microwatts  of  scattered  power  at  a  range  of  20  km.  This  power  could  easily  be  lost  in  some  atmospheric 
disturbance  or  in  transmission  or  receiving  inefficiencies.  Thus,  for  a  very  versatile  source,  the  output 
power  should  be  in  the  MW  range.  Also,  the  proposed  source  lacks  one  capability  that  would  be  necessary 
for  a  single  source  instrument  to  able  to  interrogate  the  plume  of  a  solid  rocket  motor.  That  capability  is  the 
ability  to  frequency  sweep  across  the  terahertz  regime.  Frequency  sweeping  is  needed  to  solve  the  single 
Rayleigh  theory  equation  with  the  two  unknowns  of  particle  radius  and  number  density  as  described  in  Ch. 
5  Sec.  C.  Without  this  capability,  several  sources  would  be  needed  to  conduct  an  accurate  diagnostic.  The 
source  is  only  one  issue  for  the  instrument;  an  antenna  would  still  need  to  be  developed  in  order  to  make 
use  of  the  source.  At  a  range  of  20  km  the  antenna  could  be  3  m  in  diameter.  While  not  the  largest  antenna 
by  any  measure,  it  is  still  a  sizeable  undertaking.  The  antenna  setup  would  also  need  the  ability  to  track  the 
trajectory  of  the  rocket  in  order  to  provide  continuous  data  on  the  rocket  plume  characteristics.  Also,  the 
design  for  the  notional  terahertz  instrument  relied  heavily  on  a  high  gain  to  be  realized  by  the  antenna.  The 
transmission  efficiencies  of  the  antenna  and  the  feed  antenna  would  need  to  be  adequate  in  order  to 
generate  the  necessary  gain. 

Another  problem  with  using  an  actual  terahertz  instrument  would  be  the  atmospheric  conditions 
on  the  day  of  use.  Terahertz  radiation  is  easily  absorbed  by  water.  Thus,  an  excess  of  moisture  would 
largely  hinder  the  operation  of  the  terahertz  instalment.  Figure  40  below  shows  the  attenuation  effects  of 
several  different  weather  conditions.44 


Figure  40:  Attenuation  due  to  water  for  terahertz  instrument44 
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It  is  important  to  note  that  the  y-axis  is  attenuation  per  km.  This  means  that  on  a  foggy  day  the  one-way 
loss  for  the  terahertz  instrument  would  be  about  20  dB.  This  would  absolutely  destroy  any  chance  of 
performing  any  accurate  measurements. 

The  last  major  problem  that  would  have  to  be  addressed  prior  to  conducting  actual  real-time 
diagnostics  is  the  inter-particle  interactions.  During  the  tests  of  the  proof-of-concept  experiment,  there 
were  instances  when  adding  more  scattering  centers  produced  the  appropriate  increase  in  measured 
scattered  power.  The  following  figures  are  from  an  earlier  setup  of  the  experiment  in  which  the  receiving 
horn  only  had  a  10  dB  gain  instead  of  the  14.7  dB  gain  horn  used  for  the  other  runs.  Figure  41,  Figure  42, 
and  Figure  43  show  the  effects  of  adding  more  scattering  centers.  The  spheres  used  were  the  9.53  mm  and 
4.76  mm  radii  aluminum  spheres.  The  triple  sphere  setup  is  similar  to  the  setup  seen  in  Figure  16.  The 
three  spheres  were  placed  approximately  14  cm  apart.  The  spheres  were  in  a  line  perpendicular  to  the 
direct  line  of  sight  of  the  transmitting  antenna  as  in  the  two  sphere  setup. 


—  9.53  mm 
■  4.76  mm 


Figure  41:  Scattered  power  versus  source  power  for  one  aluminum  sphere 


Source  Power  (W) 


■  9.53  mm 

■  4.76  mm 


Figure  42:  Scattered  power  versus  source  power  for  two  aluminum  spheres 
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■  9.53  mm 

■  4.76  mm 


Figure  43:  Scattered  power  versus  source  power  for  three  aluminum  spheres 

As  can  be  seen  from  the  above  figures,  the  expected  doubling  occurs  between  the  single  sphere 
and  the  double  sphere  tests.  However,  once  a  third  sphere  is  added,  the  scattered  signal  actually  decreases. 
This  seems  to  be  the  result  of  inter-particle  interactions  which  lead  to  a  lower  returned  scattered  signal. 
Drolen  and  Tien7  investigated  this  phenomenon  in  regards  to  independent  and  dependent  scattering  in 
packed-sphere  systems.  They  show  that  a  dependent  scattering  system  has  a  decreased  scattering 
efficiency.  This  may  be  the  cause  for  the  decrease  in  the  triple  sphere  scattered  power.  However,  more 
tests  would  need  to  be  conducted  in  order  to  confirm  this.  Also,  it  would  also  have  to  be  determined  if  the 
environment  inside  in  the  solid  rocket  exhaust  plume  falls  into  either  the  independent  regime  or  the 
dependent  regime  for  scattering. 
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CHAPTER  9 


CONCLUSION 

This  study  has  developed  the  notional  plans  for  a  real-time,  non-intrusive  diagnostic  tool  for  the 
characterization  of  a  solid  rocket  exhaust  plume  using  Rayleigh  scattering  theory  and  terahertz  radiation. 
This  study  has  shown  the  feasibility  of  measuring  particle  size  based  on  the  collected  Rayleigh  scattered 
light  if  the  particle  number  density,  index  of  refraction,  and  scattering  efficiency  is  known.  Terahertz 
radiation  has  captivated  many  researchers  because  of  its  relatively  unexplored  nature  and  wide  range  of 
applications.  It  occupies  the  part  of  the  spectrum  between  100  GHz  and  10  THz.  Terahertz  sources 
currently  have  relatively  low  output  powers  with  a  maximum  in  the  tens  of  Watts.  Despite  these  low 
powers,  some  terahertz  sources  have  attained  high  intensity  ratios  based  on  Rayleigh  scattering  theory.  A 
high  output  power  is  desired  for  the  projections  of  this  study.  Thus,  UCCS  researchers  are  developing  a 
compact  gyration  source  that  will  operate  at  300  GHz  with  an  output  power  of  1  kW.  This  source  can 
provide  the  necessary  output  for  the  notional  terahertz  diagnostic  instrument. 

The  optical  properties  of  alumina  were  researched  with  emphasis  put  on  the  index  of  refraction.  It 
was  determined  that  a  dependence  exists  between  the  temperature  of  the  particles  and  the  wavelength  of 
light  for  this  property.  Therefore,  studies  must  be  conducted  in  order  to  determine  the  index  of  refraction 
of  alumina  at  terahertz  frequencies  at  the  temperatures  experienced  inside  the  rocket  plume. 

The  equations  for  two  elastic  scattering  theories  were  analyzed  for  application  to  this  study.  Mie 
theory  was  briefly  discussed  and  can  be  applied  to  any  spherical  particle  with  any  incident  wavelength  of 
light.  Rayleigh  theory  which  is  an  approximation  of  Mie  theory  can  only  be  applied  if  the  particle  radius  is 
ten  times  smaller  than  the  wavelength.  The  intensity  ratios  for  both  theories  depend  on  the  particle  radius, 
the  number  density,  the  wavelength  of  the  incident  light,  and  the  index  of  refraction  of  the  material.  The 
intensity  ratio  from  Rayleigh  theory  was  resolved  into  an  equation  for  determining  the  power  received 
based  on  the  dimensions  of  the  detector  and  the  experimental  setup  along  with  the  particle  radius,  index  of 
refraction,  and  wavelength  of  light.  The  issue  of  using  one  equation  with  two  unknowns  was  discussed  and 
a  solution  could  lie  in  the  use  of  a  frequency  sweeping  source.  By  increasing  the  frequency  one  could 
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effectively  turn  off  the  scattered  power  of  larger  particles  because  of  the  nature  of  Rayleigh  and  Mie 
scattering.  Using  this  technique,  the  scattering  for  groups  of  particles  could  be  measured  and  their  average 
sizes  and  number  densities  determined. 

It  was  found  that  the  trends  of  the  data  matched  the  expected  general  trends  developed  in  the 
predictions.  Best  fit  predicted  lines  were  generated  using  Eq.  (29)  with  a  scattering  efficiency  of  100%  and 
varying  indices  of  refractions  for  the  materials.  The  indices  of  refraction  for  aluminum,  brass,  and  stainless 
steel  were  3.4,  6.0,  and  1.85  respectively.  The  percent  differences  were  reported  and  found  to  be  small  for 
the  assumed  parameters,  but  grew  with  decreasing  sphere  radius.  Many  possible  sources  of  error  were 
presented  which  may  have  contributed  to  the  differences.  The  most  prominent  of  these  error  sources  is  the 
low  signal  to  noise  ratio  produced  by  the  tests.  Also,  the  uncertainty  surrounding  the  index  of  refraction 
and  scattering  efficiency  contribute  to  the  error. 

The  proof-of-concept  experimental  results  paved  the  way  for  the  development  of  a  notional 
terahertz  diagnostic  instrument.  First,  Rayleigh  theory  was  combined  with  free-space  transmission 
equations  and  antenna  theory  for  a  parabolic  reflector  antenna.  A  general  equation  was  generated  for  the 
power  ratio  for  alumina  particles  inside  a  solid  rocket  exhaust  plume.  From  this  general  equation, 
predictions  were  made  for  the  performance  of  a  notional  terahertz  instrument.  The  importance  of 
frequency  matching  to  particle  size  was  highlighted  in  order  to  preserve  the  transmitted  to  received  power 
ratio.  Also,  the  predictions  were  shown  for  a  system  operating  at  20  km  away  from  the  rocket  plume  with  a 
3  meter  diameter  parabolic  dish  antenna.  Based  on  these  predictions  and  the  development  of  the  1  kW 
gyrotron  source,  the  notional  system  was  determined  to  be  able  to  generate  micro-Watts  of  scattered  power 
at  the  given  range  using  3  meter  diameter  transmitting  and  receiving  dishes. 

Some  complications  that  would  arise  in  the  development  and  application  of  an  actual  terahertz 
diagnostic  instrument  were  analyzed.  The  largest  complication  to  the  operation  of  such  a  system  would  be 
atmospheric  disturbances  and  attenuations.  Water  was  recognized  as  the  biggest  hindrance  to  operation. 
From  these  problems  it  is  recommended  that  an  even  higher  power  source  be  used  preferably  something  in 
the  MW  range.  Another  problem  would  be  in  development  of  the  antenna  as  far  as  construction,  tracking 
capability,  and  efficiency.  The  last  problem  encountered  is  the  effects  of  inter-particle  interactions.  Some 


Distribution  A:  Approved  for  public  release.  Distribution  unlimited. 


59 


experimental  data  was  presented  showing  the  effect  it  can  have  in  lowering  the  scattering  efficiency.  This 
issue  would  have  to  be  addressed  and  characterized  in  order  to  use  the  terahertz  instrument  effectively. 

This  study  has  developed  a  notional  non-intrusive,  real-time  diagnostic  instrument  using  Rayleigh 
scattering  and  terahertz  radiation.  It  has  also  shown  the  feasibility  of  such  a  device  through  a  proof-of- 
concept  experiment  conducted  at  microwave  frequencies  with  appropriately  scaled  particles.  While  issues 
will  need  to  be  addressed  by  future  work  in  order  to  create  and  operate  an  actual  diagnostic  tool,  many 
potential  solutions  are  presented.  A  compact  diagnostic  instrument  of  this  nature  can  have  enormous 
potential  in  the  future  of  characterizing  the  properties  of  solid  rocket  plumes. 
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APPENDIX  A:  ROCKET  PLUME  CONSTITUENT  DATA 
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APPENDIX  B:  MATLAB  CODE  FOR  EXPERIMENTAL  PREDICTIONS 


9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9-9'9'9'&-9'9'9'9'9'9'9'9-9'9'9'9' 

ooooooooooooooooooooooooooooooooooooooooooooooooo 

%  Microwave  Rayleigh  Predictions  2kW 

o, 

o 

%  This  script  calculates  the  received  Rayleigh  scattered  power  from  a 
%  variety  of  different  ball  sizes  and  materials  based  on  user  inputs. 
%  It  then  plots  the  results  as  scattered  power  vs  source  power 
%  and  scattered 

%  power  versus  particle  radius. 

o, 

o 

%  Author:  Anthony  O'Shea 

o, 

o 

%  Inputs: 

%  Parameters  of  experiment 

o, 

o 

%  Outputs: 

%  Pr  -  the  scattered  power  received 

g, 

o 

%  Globals:  None 

o, 

o 

%  Constants:  None 

o, 

o 

%  Coupling:  None 

o, 

o 

%  Refernces:  40  and  41 

9'9'9'9'9-9'9'9'9'9'9'9'9'9'9-9'S-9'9'9'S-9'9'9'9'9'9'9'9-9'9'9'9'9'9'9'9-9'9'9'9'9'9'9'9'9'9-9'9' 

ooooooooooooooooooooooooooooooooooooooooooooooooo 


%clear  all 
%clc 

%Define  variables  for  integration 
syms  theta  phi 
%Define  Sphere  Radii 

r  =  [0.75/2*2.54/100,0.625/2*2.54/100,0.5/2*2.54/100,0.375/2*2.54/100] 
%r  =  [0.75/2*2.54/100,0.375/2*2.54/100]; 

%Set  Input  Power 

Pt  =  [100,400,700,1000,1300,1600,1900,2100]; 

%Input  Experimental  Setup  Parameters 
H  =  4.875; 

W  =  6.25; 

D  =  56; 

T  =  121.75; 

L  =  122; 

Cl  =  96.25; 

C2  =  121.5; 

thetac  =  acos (DA2-C1A2-C2A2+2*C1*C2) ; 

thetaa  =  asin(W/2/L); 

phia  =  asin(H/2/T); 

phil  =  pi/2-phia; 

phi2  =  pi/2+phia; 

thetal  =  thetac-thetaa; 

theta2  =  thetac+thetaa; 

%Set  Index  of  Refraction  and  Scattering  Efficiency 
%m  =  1.39; 

%m  =  2.0016; 

%m  =  2.50; 
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m  =  1 0  A  9 ; 
eta  =  .25; 

%Input  Antenna  Gains  and  Wavelength 
dBl  =  10; 
dB2  =  14.17; 
lambda  =  0.1224; 

%Define  cross-sectional  area  of  sphere 
R  =  Cl*2. 54/100; 

%Setup  iteration  loop 
Pr  =  zeros (8,1); 
count  =  0; 
for  j  =1:4 

rp  =  r  ( j ) ; 

%Calculate  sphere  volume 
V  =  4/3*pi*rpA3; 
for  n  =  1 : 8 

Ph  =  Pt (n) ; 

%Calculate  differential  scattering  cross-section 
sigss  =  9*piA2*VA2* ( (mA2-l) / (mA2+2) ) A2* . . . 

int ( int ( ( sin (phi ) ) A3 , phi , phil , phi 2 ) , theta, thetal , theta2 ) / . . . 
lambdaA4 ; 

%Calculate  the  power  received 

%Pr(n,j)  =  10A (dBl/5) * (lambda/ (4*pi*R) ) A2*sigss/ (pi*rpA2) *Ph; 

Pr (n,  j )  =  eta*2*10A (dBl/10) *10A (dB2/10) * (lambda/ (4*pi*R) ) A2* . . . 

sigss/ (pi*rpA2 ) *Ph; 
count  =  count+1 

end 

end 

%Generate  plot  for  scattered  power  versus  source  power 
plot (Pt, Pr ( : , 1) ,  '  k'  ) 
hold  on 

plot (Pt, Pr ( : , 2) , ' — k' ) 
plot ( Pt , Pr ( : , 3 ) , ' - . k ' ) 
plot (Pt, Pr ( : , 4 ) , ' : k ' ) 

legend ( ' r  =  9.53  mm ' ,  ' r  =  7.94  mm ' ,  ' r  =  6.35  mm ' ,  ' r  =  4.76  mm '  ) 

%legend ( ' r  =  9.53  mm ' , ' r  =  4.76  mm ' ) 
hold  off 
Pr 

%Generate  plot  for  scattered  power  versus  sphere  radius 

figure  (2 ) 

plot (r, Pr (8,  :  ) ) 
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APPENDIX  C:  MATLAB  CODE  FOR  THZ  INSTRUMENT  PERFORMANCE 


ooooooooooooooooooooooooooooooooooooooooooooooooo 

%  Terahertz  Rayleigh 

o, 

o 

%  This  script  calculates  the  Rayleigh  power  ratio  from  a 
%  notional  terahertz  instrument  and  alumina  particles. 

%  It  then  plots  the  results  as  scattered  power  versus  source  power 
%  and  scattered  power  versus  particle  radius. 

o, 

o 

%  Author:  Anthony  O'Shea 

o, 

o 

%  Inputs: 

%  Parameters  of  notional  instrument 

Q_ 

O 

%  Outputs: 

%  ratio  -  the  Rayleigh  power  ratio 

g, 

o 

%  Globals:  None 

Q, 

O 

%  Constants:  None 

o, 

o 

%  Coupling:  None 

O, 

o 

%  Refernces:  40  and  41 

9'2'9'9'9'2'9'9'9'2'9'9'9'9'9'5'9'9'9'5'9'2'9'2'9'2'S-9'9'2'9'2'S-9'9'2'S-9'9'2'9'2'9'2'9'2'9'9'9' 

ooooooooooooooooooooooooooooooooooooooooooooooooo 

clear  all 
clc 

%Define  variables  for  integration 
syms  theta  phi 

r  =  100e-6;  %Radius  of  particle  (m) 
lambda  =  3e8/300e9;  %Wavelength  (m) 

V  =  (4/3) *pi*rA3;  %Volume  of  particle  (mA3) 
m  =  1.76;  %Index  of  refraction 
%Set  up  for  loop 
ratio  =  zeros ( 100 , 1 )  ; 

R  =  zeros  (100, 1)  ; 
count  =  0; 

%Define  radius  of  dish  (m) 
rdish  =  [0.5,0.75,1,1.25,1.5]; 

%rdish  =  [1.5,1.25,1,0.75,0.5]; 

Rmax  =  30000; 

%Run  iterative  loop 
for  j  =1:5 

rd  =  rdish  ( j ) ; 
for  n  =  1:100 

count  =  count+1 

%Calculate  minimum  distance  based  on  far-field  criteria 
Rmin  =  2*(2*rd)A2/0.001; 

R(n,j)  =  Rmin+ (n-1) * (Rmax-Rmin) /100; 

%Calculate  antenna  gain 

dB  =  17 . 8+20*1 og 10 (2*rd) +20*1 og 10 (300) ; 

%Determine  integral  bounds  for  Rayleigh  theory 
phil  =  pi/2-atan (rd/R (n, j ) ) ; 
phi2  =  pi/2+atan ( rd/R (n, j ) ) ; 
thetal  =  pi/4-atan (rd/R (n, j ) ) ; 
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theta2  =  pi/4+atan (rd/R (n, j ) ) ; 

%Calculate  differential  scattering  cross-section 
sigss  =  9*piA2*VA2* ( (mA2-l) / (mA2+2) ) A2*  . . . 

int ( int ( (sin (phi ) ) A3 , phi , phil , phi 2 ) , theta, thetal , theta2 ) /lambdaA4 
%Calculate  ratio  between  scattered  power  and  source  power 
ratio  (n, j )  = 

10 A (dB/5 ) * ( lambda A2 / (4*pi*R(n,j) ) ) A2*sigss/ (pi*rA2)  ; 
end 

end 

%Plot  results 

plot (R ( : , 1 ) , ratio  {:,!),  '  -  '  ) 
hold  on 

plot (R ( : , 2 ) , ratio  (  : , 2 ) ,  '  -  .  '  ) 
plot (R ( : , 3 ) , ratio  ( : , 3 ) ,  '  :  '  ) 
plot (R ( : , 4 ) , ratio  (  : , 4 ) ,  '  -  -  '  ) 
plot (R ( : , 5 ) , ratio  ( : , 5 ) ,  ' -  '  ) 
hold  off 
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